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Abstract. In this paper, we consider an incompressible viscous flow without surface tension in a finite- 
depth domain of three dimensions, with free top boundary and fixed bottom boundary. This system 
is governed by a Naiver-Stolses equation in above moving domain and a transport equation for the top 
boundary. Traditionally, we consider this problem in Lagrangian coordinates with perturbed linear 
form. In the series of papers [1], [2] and [3], I. Tice and Y. Guo introduced a new frameworlt using 
geometric structure in Eulerian coordinates to study both local and global wellposedness of this system. 
Following this path, we extend their result in local wellposedness from small data case to arbitrary 
data case. Also, we give a simpler proof for global wellposedness in infinite domain. Other than the 
geometric energy estimates, time-dependent Galerkin method, and interpolation estimate with Riesz 
potential and minimum count, which are introduced in these papers, we utilize three new techniques: 

(1) using e-Poisson integral to construct a diffeomorphism between fixed domain and moving domain; 

(2) using bootstrapping argument to prove a comparison result for steady Navier-Stokes equation for 
arbitrary data of free surface; (3) redefining the energy and dissipation to replace the original compli- 
cated bootstrapping argument to show interpolation estimate. 

Keywords: free surface, geometric structure, interpolation 

1. Introduction 

1.1. Problem Presentation. We consider a viscous incompressible flow in the moving domain, 
(l-l-l) ^{t) = {y e S X i? I - %i,y2) < ya < ?7(2/i, 2/2, 0} 

Here we can take either E = i?^ or E = (iiT) x {L2T) for which T denotes the 1-torus and L^,L^ > 
the periodicity lengths. The lower boundary b is fixed and given satisfying 

(1.1.2) < 6(yi, 2/2) < constant and heC°°{Y.) 

When E = R^, we further require that b(yi, 1/2) bo — positive constant as -t- I2/2I — > 00 and b — b^E 
H'iT) for any s > 0. When E = (LiT) x (L2T), we just denote 60 = l/2(max{6(yi, y2)}+min{6(yi, 2/2)})- 
It is easy to see this also implies 6 — 60 G H^{Yj) for any s > since E is bounded. We denote the initial 
domain 17(0) — VLq. For each t, the flow is described by velocity and pressure {u,p) : V,{t) t-^ x R 
which satisfies the incompressible Navier-Stokes equation 

dfU + u ■ Vu + Vp = /.tAu 
V • ?i = 

{pi — fM{u))iy = grjv 
u = 



(1.1.3) 



in 


n{t) 




in 


n{t) 




on 


{ys = 


»7(yi,y2,t)} 


on 


{ys = 


-b(yi,y2)} 


on 


{ys = 


'7(yi,y2,t)} 


in 


no 




on 


E 





u{t = 0) = uo 
[ r,{t = {))=r^o 

for v the outward-pointing unit normal vector on {2/3 ~ 77}, / the 3x3 identity matrix, (Dm)^ = diUj+djUi 
the symmetric gradient of u, g the gravitational constant and fi > the viscosity. As described in [1], 
the fifth equation in (1.1.3) implies that the free surface is convected with the fluid. Note that in 
(1.1.3), we have make the shift of actual pressure p by constant atmosphere pressure Patm according to 

P ^ P + 9y3 - Patni- 

We will always assume the natural condition that there exists a positive number p such that r]o-\-b > p > 
on E, which means that the initial free surface is always strictly separated from the bottom. Also without 
loss of generality, we may assume that p — g — 1, which in fact will not infect our proof. In the following, 
we will use the term "inflnite case" when E = i?^ and "periodic case" when E = (LiT) x (^2^). 



2000 Mathematics Subject Classification. 35Q30, 35R35, 76D03, 76E17. 

1 



2 



LEI WU 



1.2. Previous Results. Historically, this problem is studied in several different settings according to 
whether we consider the viscosity and surface tension, and the different choices of domains. 

For the inviscid case, traditionally, we replace the no-slip condition u = with no penetration condition 
u • n = on Efc. It is often assumed that the initial fluid is irrotational and then this curl- free condition 
will be preserved in the later time. This allows to reformulate the problem to one only on the free 
surface. The local wellposedness in this framework was proved by Wu [12, 13] and Lannes [14]. Local 
wellposedness without irrotational assumption was proved by Zhang-Zhang [15], Christodoulou-Lindblad 
[16], Lindblad [17], Coutand-Shkoller [18] and Shatah-Zheng [19]. In the viscous case, vorticity will be 
naturally introduced at the free surface, so the surface formulation will not work. Also, only in the 
irrotational case the global wellposedness was shown as in Wu [22] and Germain-Masmoudi-Shatah [23]. 

In the viscous case without surface tension, the local wellposedness of equation (1.1.3) was proved by 
Solonnikov and Beale. Solonnikov [11] employed the framework of Holder spaces to study the problem in 
a bounded domain, all of whose boundary is free. Beale [4] utilized the framework of L'^ spaces to study 
the domain like ours. Both of them took the equation as a perturbation of the parabolic equation and 
make use of the regularity results for linear equations. Abel [20] extended the result to the framework 
of LP spaces. Also, Hataya [24] proved the global wellposedness in periodic case. 

Many authors have also considered the effect of surface tension, which can help stabilize the problem 
and gain regularity. However, most of these results are related to global wellposedness for small data, 
e.g. in Bae [21]. 

Almost all of above results are built in the Lagrangian coordinate and ignore the natural energy 
structure of the problem. In [1], [2] and [3], Y. Guo and I. Tice introduced the geometric energy 
framework and prove both the local and global wellposedness in small data. The main idea in their proof 
is that instead of considering perturbation of the parabolic equation, we restart to prove the regularity 
under time-dependent basis. In our paper, we will follow this path and employ a similar argument. 

In [1], it is always assumed that the initial data uq and 770 is sufficiently small to guarantee that: (1) 
the transform between fixed domain and moving domain is a diffeomorphism; (2) the elliptic estimate in 
linear Navier-Stokes problem achieves the best regularity; (3) in proving boundedness and contraction of 
iteration sequence, the constant is small enough to get required estimate. In Section 2 of our paper, we 
will drop this assumption and prove the local wellposedness for arbitrary initial data. Hence, we have to 
use different techniques to recover the above three results. 

Also, in proving global wellposedness and decaying property in [2], the author introduced a quite 
complicated interpolation argument to increase the decaying rate of energy with minimum count. In 
Section 3 of our paper, through redefining the energy and dissipation, we greatly simplify the proof in 
that: (1) avoiding the complicated bootstrapping argument for interpolation results; (2)avoiding the 
bootstrapping argument in comparison theorem; (3) handling the energy and dissipation 1-minimum 
count and 2-minimum simultaneously instead of separately. 

1.3. Geometric Formulation. In order to work in a fixed domain, we want to flatten the free surface 
via a coordinate transformation. Beale introduced a flatten transform in [4] and we will use a slightly 
modified version as the geometric transform. We define a fixed domain 

(1.3.1) = {a; e S X i? I - 60 < .X3 < 0} 

for which we will write the coordinate x G In this slab, we take S : {xz = 0} as the upper boundary 
and Ef, : {2:3 = —60} as the lower boundary. Simply denote x' = {xi,X2) and x = {xi,X2,X3). In order to 
introduce the mapping between fl and 0(t), we need to utilize slightly different techniques to construct 

the extension of free surface. 

1.3.1. Geometric Formulation for Local Wellposedness. We will prove the local wellposedness for arbi- 
trary initial data and arbitrary smooth bottom, so the mapping should be related to the data itself. 
Hence, we deflne the e-Poisson integral for jj: 

f /fl2 57(Oe''^'''='e2'^'''''^d^ for infinite case 

(1.3.2) fj'=W=\ 

[ Ene{L-'z)xiL-^z) e2™-^'e^l"l^=*??(n) for periodic case 

where fi{^) and r/(n) denotes the Fourier transform of ri{x') in either continuous or discrete form and 
< c < 1 is the parameter. The detailed definitions are shown in (2.1.1) and (2.1.10). 
Consider the geometric transform from f2 to ^{t): 

b x^ 

(1-3.3) : {xi,X2,X3) ^ {xi,X2, —x^+fj^l + — )) = (1/1,2/2,2/3) 
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This transform maps Cl into Sl{t) and its Jacobian matrix 



(1.3.4) 



and the transform matrix 




(1.3.5) 



where 



(1.3.6) 



1 -A^K^ 
1 -B^K" 




6=1 + 



X3 

bo 



A' = + difj'b 

b fj^ 
J' = - + J-+dsf]'b 
bo bo 



bo 



Note that this transform is not necessarily a homomorphism for arbitrary e. By our assumption on initial 
data that ??o(2;') + 60(2;') > p> and theorem 2.7, wc can always choose a sufficiently small e depending 
on ||?7o||^5/2 such that there exists a (5 > satisfying J'^(O) > (5 > 0, then this nonzero Jacobi implies the 
transform is a homomorphism at f = 0. If we further assume r\o € we can conclude is a 

diffeomorphism from <1 to J7o- In the following, we will just write fj instead of ff for simplicity, and 
the same fashion applies to A, A, B, J and K. 
Define some transformed operators as follows. 



(1.3.7) 



^Af = ■ ^Af 
J\f = (-9177,-5277,1) 
(D^M)y = AikdkUj + AjkdkUi 



where the summation should be understood in the Einstein convention. If we extend the divergence 
Vyi- to act on symmetric tensor in the natural way, then a straightforward computation reveals that 

• Sa{p, u) = ^aP — ^A'u for vector fields satisfying ■ u = 0. 
In our new coordinate, the original equation system (1.1.3) becomes 

dtu — dtfjbKdsu + u ■ W^u — A^u + V^p = in f2 

• u = in n 

SA{p,u)Af = rjN' on E 

u — on Eb 

u{x,0)=uo{x) in fl 



(1.3.8) 



dtrj + uidiTj + U2d2r] ■ 

77(2;', 0) =7,0 (.t') 



W3 



on E 
on E 



where we can split the system into a Naiver-Stokes equation and a transport equation. 
Since A is determined by r? through the transform, so all the quantities above is related to r?, i.e. the 
geometric structure of the free surface. This is the central idea of our proof. It is noticeable that in 
proving local wellposedness of above equation system, we must verify <^{t) is a diffeomorphism for 
any t € [0, T], where the theorem holds. 



1.3.2. Geometric Formulation for Global Wellposedness. In the global part, we need to consider the 
infinite case when domain Q,{t) always possess a fiat bottom, which means b = constant. Also our 
theorem only holds for sufficiently small data, so we may directly define f} as the Poisson integral of rj 
which is shown in (A.3.1), i.e. ry = Vrj. Then the map can be taken as 



(1.3.9) 



$ : {xi,X2,Xi) ^ {xi,X2,Xi+f}{l + —)) = (2/1,^2,2/3) 
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which also successfully transform ft to n{t). It is obvious that the Jacobi matrix and transform matrix 
are identical to those in the local case except that we need to redefine the quantities 

b=l + x^/b 

(1.3.10) A = difjb B = d^nb 

J =l + 7]/b + d3fib K=l/J 

Naturally, the definition of weighted operators in (1.3.7) and the transformed equation (1.3.8) are the 
same as those in local part with replacing the related quantities with above definition. We can easily 
see that as long as ||%||7j5/2(j]) is sufficiently small, $ is a diffeomorphism. Hence, we do not need the 
discussion as local part. 

1.4. Main Theorem. In this paper, we will prove both the local wellposedness and global wellposedness 
for higher order regularity. 

Theorem 1.1. Let N > 3 be an integer. Assume the initial data r]o + b > 26o^ > for some ^ > 0. 
Suppose that Uq and t]o satisfy the estimate Kq ~ |luo|l^/2jv + |1?7oI1h2«+i/2(x;) < oo as well as the N*^ 
compatible condition (2.6.3). Then there exists < Tq < 1 such that for <T < Tq, there exists a unique 
solution {u,p,ri) to system (1.3.8) on the interval [0,T] that achieves the initial data. Furthermore, the 
solution obeys the estimate 



(1.4.1) 



N 

E 

j=0 



sup 

0<t<T 



dlu 



H2N-23 



+ 



N-1 

E 

3=0 



sup 

0<t<T 



N 

E 

3=0 

dip 



dlu 



fJ2N-2j + l 



N 



JJ2N-2j-l 



+ sup ||7?||^2iv+i/2(s) +E ^t'^ 
\0<t<T ^ ' 0<t<T 



E 

3=0 



Jo 



. 0<t<T 
i-T 



il'?llj/2N + l/2(s) + 



Jo 



-i/2(E) 



N + 1 

+ E 

j=2 



0<t<T 

T 







2 


r 


dip 




Jo 




H2N-2j 




2 




91 V 






JJ2N 


-23+3/2 ( 







dlr] 



JJ2N -23+5/2 (J2) 

< C{flo,6)P{Ko 



where C(r2o,(S) is a positive constant depending on the initial domain D,q and separation quantity S and 
P(-) is a single variable polynomial satisfying P(0) = 0. The solution is unique among functions that 
achieves the initial data and for which the left hand side of the estimate is finite. Moreover, ri is such 
that the mapping ^{t) defined by (1.3.3) is a C'^^~'^ diffeomorphism for each t G [0,T]. 

Remark 1.2. In above theorem, H'^ denotes the usual Sobolev space in O and iI*'(S) denotes the usual 
Sobolev space on S, while we will state the definition of ||-||_^.. later. 

Because the map $(i) is a c2iv-2 diffgofnorphism, then we may change variable to y £ il{t) to produce 
solution of (1.1.3). 

Theorem 1.3. Let N > A be an integer. Suppose the initial data {uo,r]o) satisfy the appropriate com- 
patible conditions stated in local wellposedness theorem. Then there exists a k> such that if 



(1.4.2) 



f2iv(0)+ J"2iv(0) < K 



then there exists a unique solution {u,p,r]) on the interval [0,oo) that achieves the initial data. Also, the 
solution obeys the estimate 

(1.4.3) gjv(oo) < C(f2w(0) + J-2iv(0)) < Ck 

where C > is a universal constant. 

Remark 1.4. The detailed definition for £2n, Tin and Qn is presented in (3.2.1) to (3.2.12). This 
theorem easily implies the algebraic decaying of energy £n+2,2 and £n+2,\- 

1.5. Convention and Terminology. We now mention some of the definitions, bits of notations and 

conventions we will use throughout the paper. 

(1) We will employ Einstein summation convention to sum up repeated indices for vector and tensor 
operations. 
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(2) Throughout the paper C > will denote a constant only depend on the parameter of the prob- 
lem, N and f2, but does not depend on the data. They are referred as universal and can change 
from one inequality to another one. When we write C{z), it means a certain positive constant 

depending on quantity z. 

There are two exceptions to above rules in the local part. The first one is that in the elliptic 
estimates and Korn's inequality, there are constants depending on the initial domain CIq- Al- 
though this should be understood as depending on the initial free surface rjo, since its dependent 
relation is given implicitly and cannot be simplified further, we will also call them universal. The 
second one is that we will also call the constant C{6) universal, though 6 is determined by initial 
domain Hq. To note that apart from these, all the other constants related to initial data CIq, uq 
and rjo should be specified in detail. 

(3) We will employ notation a < 6 to denote a < Cb, where C is a universal constant as defined 
above. 

(4) We will write P{-) to denote the single variable polynomial. This type of polynomial always 
satisfies P(0) = 1 unless it is specified as an exception. This notation is used to denote some 
very complicated polynomial expressions, however whose details we do not really care. This kind 
of polynomial may change from line to line. 

(5) For convenience, we will typically write = L^, except for notation L'^{[0,T]; H''). We write 
H''{n) with fc > and iJ"(S) with s e R for standard Sobolev space. Same style of notations 
also holds for L^{[0,T]- H'') and L°°([0, T]; ff*^). When we write \\-\\Hk, this always means the 
Sobolev norm in ft, otherwise, we will point out the exact space it stands for, e.g. |H|//fc(s)- A 
similar fashion is adopted for j|-j|j^2j/fc and ||-|j^cx=jjfc. 

(6) We write the multi-indices ^1+™ = {a = (ao, ai, . . . , am)} to emphasize that the 0-index term 
is related to temporal derivatives. For just spatial derivatives, we write without the 0-index. 
We define the parabolic counting of such multi- indices as \a\ = 2aQ + ai + . . . + a^- We always 
write Df to denote the horizontal derivative of / and V/ for the full derivative. 

For a given norm || || and integers fc, m > 0, we introduce the following notation for the sums of 
spatial derivatives. 

\\Dtf\^= E iiv^/ir= E 

a&P m<\a\<k aefP m<\a\<k 

QeNi+2 m<|Q|<fe QeNi+3 'm<\a\<k 

where D and V means horizontal derivatives and D and V means full derivatives. 
Also we define 

\\D^f\^ = \\Dlff llv'^/f = ||V^/f 

ii^vir = ii5,vir iivvir=iiv^/ir 

1.6. Structure of This Paper. In Section 2, we will discuss the local behavior for arbitrary initial data, 
which includes three main parts: estimates for linear Navier-Stokes equation, estimates for transport 
equation and the construction of iteration. Combining all these, the local wcllposcdness easily follows. 
In Section 3, we will prove the global wellposedness for small data in horizontally infinite domain. Since 
many lemmas have been proved in [2], here we only present the main part of our simplification in 
interpolation theorem and comparison theorem, with merely stating the result for other parts. 



2. Local Wellposedness for Large Data 

2.1. e-Poisson IntegraL In this section, we will define the e-Poisson integral, which is utilized to 
construct geometric transform as (1.3.3), and discuss the basic properties of it. 

2.1.1. Poisson Integral in Infinite Case. For a function / defined on E = B?, the e-Poisson integral is 
defined by 

(2.1.1) V'f{x',X3)= [ /(Oe^l«l"^e2--'-«dC 

where /(^) denotes the Fourier transform of f{x') in and < e < 1 is the parameter. Although V^f 
is defined on R'^ x (— oo,0), we will only consider the part R"^ x (— 6o,0) here. 
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Lemma 2.1. Let V^f be the e-Poisson integral of function f which is in homogeneous Soholev space 
jj9-i/2(s) /or g e N. Then we have 



(2.1.2) 

where C > is a constant independent of e. In particular, we have 
(2.1.3) 

Proof. By definition of Fourier transform, Fubini theorem and Parseval identity, we may bound 



V^T'VllLo < i'^^ f I /" l?l'V"(0'e2^'«'"^d^dx3<(27r)2« / |^|'« /(O ' T ' 



IR^ J -bo 

< 7 ll/lli/<!-l/2(E) 



2 ^1 _ g-2e6o|C| 



(27r) 



2? 



2e 



i^r^"' fio 



R^ 



dC 



We can simply take C = tt to acliieve the estimate. Furthermore, considering 



1 _ e-2«bol«| 



< min 



1 



and < e < 1, we have 
(2.1.4) 

Hence, (2.1.3) easily follows. 

We also need the L°° estimate for Poisson integral. 
Lemma 2.2. Let V^f he the e-Poisson integral of function f. Then for q & N and s > 1, we have 
(2.1.5) l|V''7'^/lli^ <C||/||^,+.(5.) 

Proof. A simple application of Sobolev embedding and lemma 2.1 reveals 



□ 



□ 



The following lemma illustrate the specialty of derivative in vertical direction. 
Lemma 2.3. Let V^f be the e-Poisson integral of function f . Then we have 

(2.1.6) ||53P7lli=o <Ce 11/11 
where C > is a constant independent of e. 

Proof. We can simply bound as follows. 

(2.1.7) PsPVIlioo < C Wd^V^fWl^ < Ce' WV^fWls < Ce ||/||^5/2(s) 

Note that the first inequality is based on Sobolev embedding theorem and the third one is a straightfor- 
ward application of lemma 2.1, so we only need to verify the second inequality in detail. 
By definition, it is easy to see 



(2.1.8) 

Hence, we have 



daV'f = e 



Jr^ 



d^ 



WdsV^fWlo = ^'J^J[ I^l1/(^)r«'''^'''dedx3 = ^((27r)2|^^y° \^\'\m 



d^dxj 



< 



(2^)2 

Similarly, we can show that for i,j = 1, 2, 3 

WdsiV'fWlo < 



II^ii^'VIIho 
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So we have 
(2.1.9) 



□ 



2.1.2. Poisson Integral in Periodic Case. Suppose S = (iiT) x (L2T), we define the e-Poisson integral 
as follows. 



(2.1.10) 

where 
(2.1.11) 



E 



^2Trin-x' e\n\x3 



fin) 



ne{L-^I,)x{L-^2 



f{n)= [ f{x'y—— 



-2TTin-x' 



-dx' 



Lemma 2.4. Let V^f be the e-Poisson integral of function f which is in homogeneous Sobolev space 
ij-9-V2(5]) for qeN. Then we have 



(2.1.12) 

where C > is a constant independent of e. In particular, we have 
(2.1.13) 



Proof. By Fubini theorem and Parseval identity, we may bound 

/■O 

|2 



IIVPVIIho < (27r)2« 



S / 



^|2g ^2e\n\x3 



{27r 



|2g 



E 



n 



bo 



2q 



fin) 



dx3 



fin) 



2 / 1 — e"^'^^"!' 



2e|n| 



< 



(27r) 



29 



2e 



E 



n 



2g-l 



/(") 



< 



^ M j-m2 



lff9-l/2 



We can simply take C = tt to achieve the estimate. Furthermore, considering 



I _ g-2e6o|"| 

2e|n| 



< min 



1 



2e\n\ 



\\r^f Wlo < ^ ||/||^o(s) 



and < e < 1, we have 
(2.1.14) 

Hence, (2.1.13) easily follows. 

Similar to infinite case, we give an L°° estimate for Poisson integral. 
Lemma 2.5. Let V^f be the e-Poisson integral of function f. Then for g e N and s > 1, we have 

||v«P7llioo<c||/||UME) 

Proof. A simple application of Sobolev embedding and lemma (2.4) reveals 

(2.1.15) W^'V^fWloo < C II V«PVIIh=+V2 < C ||/||U3(E) 

We still need the following lemma to illustrate the specialty of derivative in vertical direction. 
Lemma 2.6. Let V^f be the e-Poisson integral of function f. Then we have 

(2.1.16) WdsV'fWloo <Ce\\f\\ 
where C > is a constant independent of e. 



□ 



□ 
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Proof. We can simply bound as follows. 

(2.1.17) WdsV^fWloo < C WdsV'fWl. < Ce^ WV'fWls < Ce \\f\\ 

Note that the first inequality is based on Sobolev embedding theorem and the third one is a straightfor- 
ward application of lemma (2.1), so we only need to verify the second inequality in details. 
By definition, it is easy to see 



(2.1.18) 



Hence, we have 



E 



bo 



fin) 



bo 



fin) 



< 



{2n) 



Similarly, we can show that for i,j = 1, 2, 3 



So we have 
(2.1.19) 



\\dsijV^f\\lo<-^\\dniV^ffHO 



□ 



2.1.3. Homomorphism of the Geometric Transform. 

Theorem 2.7. // the initial data satisfies 770 + &o > P > for all x' G S, then we can choose a e > 
such that there exists a 6 > satisfying J^{0) > 6 > 0. 

Proof. Naturally there always exists 6 > 0, such that r]o{x') + bo{x') > 2bo5 > for arbitrary x' G S. 
Based on lemma (2.3) and lemma (2.6), if we take e < (5^/(4C^ ll%llff5/2), we have 



(2.1.20) 



< S/2 



where ijq denotes r]^att = 0. Furthermore, by fundamental theorem of calculus, we can estimate for any 
2:3 e [-^0,0] and a;' e S 



(2.1.21) 

Since f]Q{x',0) = r]o{x'), we have 
J^(0) = 



f}'oix',X3)-fj'oi^',0)\< I \d3f]'oix',z)\dz<bo\\d3f]'o\\L^<boS/2 

bo 



(2.1.22) 



> 25 - (5/2 - (5/2 = (5 > 



This means for given initial data if we take e = S /{AC ||77o||jj5/2), then J"^ > S > 



□ 



2.2. Preliminaries. 
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2.2.1. Transform, Estimates. In order to study the linear problem in the slab domain, we will employ 
the idea that transforming the variable-coefhcient problem into a constant-coefficient problem through 
diffeomorphism. So before estimating, we need to confirm the mapping $ is an isomorphism from O to 
Q.{t) and determine the relation of corresponding norms between these two spaces. 

Lemma 2.8. Let ^ : ^ Q.' he a diffeomorphism satisfying e Hj"^^ , the Jacobi J = det(V4') > 
S > a.e. in n and - I £ H''{n) for an integer k > 3. If v £ H"^{n'), then v o & H"'{n) for 
m = 0, 1, . . . , A; + 1, and 

(2.2.1) \\v o < P{\\V^ - I\\HHn))MH-in') 

for a polynomial P{-). Similarly, for u G iJ™(0), u o g H"^{V,') for m = 0, 1, . . . , fc + 1 and 

(2.2.2) \\u o *-1||h".(0') < P(||V* - 

Let E' = *(E) denote the upper boundary of n' . If v € if™-i/2(i]') for m = l,...,k - 1, then 
vo^G if™-V2(s), and 

(2.2.3) 11^; o *||h--i/2(s) < C(||V* - I\\H^a))\MHm-U2^j:') 
//ueif'"-i/2(S)/orm = l,...,fc-l, i/ienuo*-i e ^'"-^/^(E') and 

(2.2.4) ||uo 'J'-i||j,„_i/2(s') < C(||V* - I\\HHQ))\MHm-i/2^j:) 
where C(||V^ — -'^11^^(0)) denote a constant depending on — I\\H''{n)- 

Proof. The proof is almost identical to that of lemma 2.3.1 in [1]. It is easy to sec from the original 
proof, in (2.2.1) and (2.2.2), the constant C is actually a polynomial of HV^* — IWHi'iQ)- other 
part of the proof is exactly the same. □ 

Remark 2.9. Based on our assumptions on fj and b, it is easy to see $ defined in 1.3.3 is a 
diffeomorphism satisfying the hypothesis in lemma 2.8. Moreover, the universal constant C now can be 
taken in a succinct form 

(2.2.5) C = C(||77||^.+V2(s)) 

// we only consider the estimate in the domain f2 and need to specify the constant, we can take the explicit 
form 

(2.2.6) C = P(||r?||^.+i/2(s)) 

2.2.2. Functional Spaces. Now we introduce several functional spaces. We write H''{Q) and il'^(E) for 
usual Sobolev space of either scalar or vector functions. Define 

W{t) = {u{t) e H'in) : «(t)|s, = 0} 
^ ' X{t)^{u{t)eH\n): =0,V^(i).^i(t) = 0} 

Moreover, let L'^{[0,T]; H''{n)) and L'^{[0,T]; H''{J:)) denotes the usual time-involved Sobolev space. 
Define 

o^ w = {ue L^io, nn^n)) = 0} 

^ ' X = {ueLH[0,T]:H\n)): = 0, • « = 0} 

It is easy to see W, X, W and X are all Hilbert spaces with the inner product defined exactly the same 
as in i7^(r2) or L^{[0, T]: (il)) . Hence, in the subsequence, we will take the norms of these spaces the 

same as || • or || • \\l^h^- 

Furthermore, we define the functional space for zero upper boundary. 

(2.2.9) V{t) = {u{t) e H\n) : w(t)|s = 0} 

Also we will need an orthogonal decomposition i/"(il) = Y{t) Y-^{t), where 

(2.2.10) Y^{t) = {V^(,)^(t) : ^{t) e V{t)} 

In our use of these spaces, we will often drop the {t) when there is no potential of confusion. 
Finally we will define the regular divergence- free space in H^{fl). 

Xo = {ueH\n): =o,v.« = o} 
^ ' X„^{ue L2([o, T];H\n)) : u\^, = 0, V • u = 0} 

Sometimes, we will consider the functional space with the similar properties as defined above in trans- 
formed space SI'. If that is the case, we will employ the notation W{fl') or X{Sl') to specify the space 
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they live in. The same fashion can be apphed to other hmctional spaces. 
We will use (•, ■)na to denote the normalized inner product in i.e. 

(2.2.12) {u,v)uo = / Ju-v 

Jn 

and (•, •)£2^o to denote the normalized inner product in L'^{[0,T]-, H^), i.e. 

(2.2.13) {u,v)c2no= [ [ Ju-v 

Jo Jn 

Also we will use (•, ■)w and (•, •)w to denote the regular inner product in if both arguments are in 
W or W. Moreover, we define the regular inner product on S as follows. 

(2.2.14) {u,v)-uo(i:) = j^u-v 

(2.2.15) {u,v)c-^u°(T.)= / ju-v 

Jo Jt, 

The following result gives a version of Korn's type inequality for initial domain f^o- Here we employ 
Beale's idea in [4]. 

Lemma 2.10. Suppose Oq is the initial domain and r]o € iI^/^(S). Then 
(2-2.16) Ml^no) ^ l|ID'^ll?fO(Oo) 

for any v G H^{flo) o,nd v = on {ys = —b}. 

Proof. In the periodic case, S is bounded and v\y:^ = 0, so naturally Korn's inequality is valid (see proof 
of lemma 2.7 in [4]). Then we consider the infinite case. First we prove the decaying property of initial 
surface rjo, i.e. for \a\ < 1, — >■ as |x| ^ oo in the slab. For horizontal derivative 9", 

a«?7o = / (27riO"e^^^l«le2"^'-«%(OdC 
Jr^ 

where rfo is the Fourier transform of r/o in R'^. Then 

/ (27riO"e-3l«l,,o(^) < / |Cre-3l?l|^^(^)|(i + |^|5/4^/' ^-—)d^ 

Jb? Jr^ Vi + |^| ' / 

1/2 / r e2^^3|«ld^ \ 1/2 



< 



/^jcri%(oi'(i + i^i'/*m 



The last inequality is valid since 0:3 < 0. Hence, (27ri^)"e'^'^3l^l?7'o(0 G L^{B?). A similar proof can justify 
the result if d°' = 83. Since |a;| — )• 00 naturally implies \x'\ — )• 00 in the slab, the Riemann-Lebesgue 

lemma implies |(?"?7o| — > as \x\ — > 00. 

We construct a mapping a = $(0) as defined in (1.3.3), which maps Q = {R^ : —bo < 2:3 < 0} to Qq- 
Denote a{x) = x + a{x). The above decaying property and our assumption on b leads to d"a{x) ^ as 

|a;| — >■ 00 for |q:| < 1. 

We partition the slab O into cubes 

Qj,k = {x : j < xi < j + 1, k < X2 < k + 1, -60 < X3 <0} 

In each cube, we have Korn's inequality 

l|w||^i(Q) < c(||iD)w||?jo(Q) + ||w|Iho(q)) 

Employing the compactness argument as Beale did, under the condition w = on {xs = —60}, we can 
strengthen this result to 

II^^IIhi(q) < C'IIid'^^IIho(q) 

This argument relies on the fact that for such v 

II^IIhO(Q) < C||Pt;||^0(Q) 
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Now suppose _D C ilo is the image of Tinion of such cubes contained in {|a;| > R} for some R. Applying 
last estimate to v = uo a and then transform to Oq, we have 



u 



m(D) - ^ \\^''^\\h'>{d) + ^ 



where e < 1 when R is large enough based on the decaying property. Then we use Korn's inequality in 
the bounded domain Qq — D and combine with above estimate to obtain the estimate required. □ 

Next, we will show the equivalence of certain quantities. Although the result is quite obvious for 
small initial data (see lemma 2.1 of [1]), now we need to utilize above Korn's inequality to show this for 

arbitrary data. 

Lemma 2.11. There exists a < fo < 1; such that if \\rj — 77o||j:^5/2 < eo, then the following relation 
(2.2.17) < ^ J \U\' < (1 + ||%|lHa/.(s)) ll^ll^o 

(2-2.18) I ^ Ml, < [ J \B^uf < (1 + ||r7o||«5/.(s))' IkllL^ 

holds for all u G W. 

Proof (2.2.17) is just a natural corollary of relation S < \\J\\loo < (1 + ||V^||^oc) < (1 + ||»7|lij6/2) < 
(1 + ||?7o||i/5/2) based on our assumption and Sobolev embedding. 
To derive (2.2.18), notice that 

(2.2.19) / J|D^m|^= / J|B^„u|V / J(e^M + D^„u) : (D^u - D^^m) 

Jn Jn Jn 

For the first term on the right-handed side, based on the integral substitution and Korn's inequality we 

proved above, we have 

/ J|D^„n|'dx > — ^-i /■ Jo|B^^^|2dx=-— -i [ \Bvfdy 

^ 1 II ||2 1 II ||2 

For the second term on the right-handed side, using Korn's inequality in slab ft, we can naturally estimate 

[ J{B^u + B^,u):iDAU-Dy^„u) < \\J\\^^ \\A + Aoh^ \\A- Aoh^ [ l^uf 
Jn Jn 

For given initial data, we can always take eo sufficiently small to absorb the second term into the first 
one. Then we have 

jmAu\^> ^ 



In U + ll'?o||i/5/2(s)j 

The first part of (2.2.18) is verified. For the second part, notice that 



J \Dj^u\' < (1 + ||r?||^./2) / \B^u\' < (1 + ||,yo||^a/2) max{l, \\AK\\l^ , \\BK\\l^ , 

Jn 



n 



Since it is easy to see 

max{l, \\AK\\l^ , \\BK\\l^ , ||i^||^^} < 1 + (1 + ||Vr?||^^) < (1 + ||%||hV2(e))' 

then the second part of (2.2.18) follows. 



□ 



Remeirk 2.12. Throughout this section, we can always assume the restriction of rj depending on eo is 
justified, and finally we will verify this condition for t G [0, T] in the nonlinear part. 

Now we present a lemma about the differentiability of norms in time-dependent space. 

Lemma 2.13. Suppose that u G W and dtu G W*. Then the mapping t — )• ||w(i)||jj-o is absolute 
continuous, and 

(2.2.20) ^ 11^(4)11^, = 2{dtu{t), u{t)) + \u{t)\' dtJ{t) 
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Moreover, u G C°([0,T]; iJO(O)). Ifv€W and dfV e W* as well, we have 
(2.2.21) / u{t) ■ vit) = {dtu{t), vit)) + {dtvit),u{t)) + I u{t) ■ v{t)dtJ{t) 



Proof. This is exactly the same result as lemma 2.4 in [1], so we omit the proof here. □ 
Next we show the estimate for H~^/'^ boundary functions. 

Lemma 2.14. If v e H° and W a ■ v e H° , then v ■ Af e iI-V2(S), v ■ u e H-'^/'^{T,b) (v is the outer 
normal vector on Sb ) and satisfies the estimate 

(2.2.22) \\v ■ A^||h-i/2(e) + \\v ■ z^IIh-i/^(e,) ^ (1 + MH^/^f ( Mho + \\Va ■ v\\ho ) 

Proof. We only need to prove the result on S; the result on S;, can be derived in s similar fashion 

considering JAe^ = v on Hi,. 

Let (j) e iJ^/^(5]) be a scalar function and let ^ e F be a bounded extension. We have 

/ 4)v-J\f = / JAijVi^{ej -63)= / • {v4>)J = / ^Va -vJ + v- Va^J 

< (1 + hWllHV.) ( ll'^ll^^ l|Vx • ^;||^o + Mho 

<{^ + \m\\H^f^nMHo + \\'^A-v\\Ho)\m 

Since ^ is arbitrary, our result naturally follows. □ 

Remark 2.15. Recall the space y{t) C . It can be shown that if v € 3^(i); then Va ■ v — in the 
weak sense, such that lem,m,a 2.14 implies that v-M G iJ~-^/2(S) andv-u € H~^/'^{Yji,). Moreover, since 
the elements ofy{i) are orthogonal to each for (j) G V{t), we find that v ■ v = Q on S;,. 

We want to connect the divergence-free space and divergence- ^-free space, so we define 

/ K 

(2.2.23) M = M{t) = KV^ = ^ 

\ AK BK 1 

Note that M is invertible and = JA^. Since dj{JAij) = for i = 1,2,3, we have the following 

relation 

p = • w ^ = JVa ■ V = JAijdjV, = dj{JAtjVi) ^ dj{JA^v)j = dj{M^^v)j = V • (M^^v) 

Then if • u = 0, we have V • {M~^v) = 0. M induces a linear operator Mt : u — )■ Mt{u) = M{t)u. 
It has the following property. 



Lemma 2.16. For each t G [0, T], k = 0,1,2, Mt is a bounded linear isomorphism from il^(O) to 

2(S)) 



II^{il) and from Xq to X . The bounding constant is given by (1 -|- ||??(i)||j:/9/2(2))^- U ^6 further define 



M defined by Mu{t) = Mtu{t), then it is a bounded linear isomorphism from L2([0, T]; i/'^(0)) to 
i^([0, T]; i?'^(ri)) and from Xq to X . The bounding constant is given by {1 + supQ^f^j' \\ri{t)\\ H9/2(^^-j)'^ . 

Proof. It is easy to see for each t e [0, T], 

WMtuWn. < \\Mt\\c2 \\u\U < (1 + ll^(i)ll^^./.(E))' 

for k = 0,1,2, which implies Mt is a bounded operator from H'' to H''. Since Mt is invertible, we 
can estimate ||A^~-^u|| ^ (1 + ll'7(^)ll/f9/2(E))^ Iklli/fc- Hence Mt is an isomorphism between H'^. Also 
above analysis implies Mi maps divergence- free function to divergence- free function. Then it is also 
an isomorphism from Xq to X{t). 

A similar argument can justiiy the case for -£/^([0, T]; TJ*^) and X. □ 

2.2.3. Pressure as a Lagrangian Multiplier. It is well-known that in usual Navier-Stokes equation, pres- 
sure can be taken as a Lagrangian multiplier. In our new settings now, proposition 2.2.9 in [1] gives 
construction of pressure from transformed equation 1.3.8, which is valid for small free surface. So our 
result in arbitrary initial data is not completely ready, which we will present here. 
For p e H^, we define the functional St G W* by St{v) = {p,Va ■ v)-ho. By the Riesz representation 
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theorem, there exists a unique Qt{p) E W such that St{v) = {Qt{p),v)w for all v S W. This defines a 
linear operator Qt : H'^ — > W, which is bounded since we may take v = Qtip) to see 

\\Qt{p)\\w = St{v) - {p,Va ■ v)no < (1 + l|rKi)ll^f5/2(s)) Ibll^o ||V^ • v\\^„ 

< (1 + Mt)\\HW2^s)f \\p\\m Mw = (1 + hWllH5/2(s))' WpWho \\QiXp)\\w 
so we have ||Qt(p)||i4/ ^ (1 + ll'7(^)ll^fs/2(5]-))'' IIpIIh"- Similarly, for S G W*, we may also define a bounded 
linear operator Q : L'^{[0,T]; H^) — )• W via the relation (p, • v)c2-h° = {Q{p)jv)w = ^{v) for all 
veW. Similar argument shows ||Q(p)||w < (1 + supo<t<T \\v{t)\\m/2(^sjf Wph^m- 

Lemma 2.17. Let p G H^, then there exists a v G W such that ■ v = p and \\v\\^r ^ (1 + 
11^(0 II /f.V2(s))^ Ibllff" • V instead p € L'^{[0,T]; H°), then there exists a v gW such that ■ v = p for 
a.e. t and ||v||w ^ (1 + supo<t<T ll»?Wllij5/2(s))'' ibllL^jjo- 

Proof. In the proof of lemma 3.3 of [4], it is established that for any q E the problem V • m = g 

admits a solution u G W such that ||u||j:^i < ||9||//o- A simple modification of this proof in infinite case 
can be applied to periodic case. Let define q = Jp, then 



\l\\lo= I M'= / bl'^'<ll^lli»|W|?,„<(l + ||r,(t)||^./2(s))^|b| 



HO 



Hence, we know v — M{t)u e W satisfies V^w = p and 

|2 ^ /I , ii_/^Mi n4ii_.ii2 



Mw < (1 + hWllH«/2(E))' \H\w ^ (1 + lbWlli/v.(E)r Ikll^o 



< 



(l+ll'?(*)llffO/2(s))'^lb|lHO 



A similar argument may justify the case for W. □ 

With this lemma in hand, we can show the range of operator Qt and Q is closed in W and W. 
Lemma 2.18. R{Qt) is closed in W, and R{Q) is closed in W. 

Proof For p G H'^ , Yet v gW be the solution of ■ v = p provided by lemma 2.17. Then 

ibilno < {P, • v)no = (Qt(p), v)w < \\Qtip)\\w Mw ^ \\Qtip)\\w (1 + imWHV^j:))' Mm 
such that 

^ nJn ^ \mp)\\w ^ Mho < (1 + IbWII^V^tE))' \\Qt{p)\\w 

Hence R{Qt) is closed in W. A similar analysis shows that R{Q) is closed in W. □ 

Now we can perform a decomposition of W and W. 

Lemma 2.19. We have that W = X®R{Qt), i.e. = R{Qt). Also, W = X®R{Q), i.e. = R{Q). 

Proof. By lemma 2.18, R{Qt) is closed subspace of W, so it suffices to show R{Qt)'^ = X. Let v G 
R{Qt)^, then for all p G we have 

{p,Va ■ v)uo = {Qtip),v)w = 

and hence V,a ■ v = 0, which implies R{Qt)'^ C X. On the other hand, suppose v G X, then • u = 
implies 

{Qt{p),v)w = {P^A ■ v)uo = 

for all p G . Hence V G R{Qt)'^ and we see X C R{Qf)-^. So we finish the proof. A similar argument 
can show W = A"® i?(Q). □ 

Proposition 2.20. If X G W* such that X{v) = for all v G X, then there exists a unique p{t) G 
such that 

(2.2.24) b(i), • v)ho = X{v) for all v gW 
and IbWIIno < (1 + lb(*)llH./.(E))'l|A|k- 

If A G W* such that A{v) = for all v G X, then there exists a unique p G ^^([0, T]; H^) such that 

(2.2.25) {p, SI A ■ v)c2n° = A(t^) for all v gW 
and \\p\\l2ho ^ (1 + supo<t<T lb(i)llH9/2(s))^||A||w 
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(2.3.1) 



Proof. If \(v) = for all v G X , then the Ricsz representation theorem yields the existence of a unique 
u e such that X{v) — {u,v)w for all v G W. By lemma 2.19, u = Q{p) for some p G Then 
A(w) = (Qt(p), t') w = {p{t), ■ v)w for all v&W. 

As for the estimate, by lemma 2.17, we may find v & W such that ■ v = p and ||u||^ ^ (1 + 
IIJ^IIh"- Hence, 

So the desired estimate holds and a similar argument can show the result for A. □ 

2.3. Elliptic Estimates. In this section, we will study two types of elliptic problems, which will be 
employed in deriving the linear estimates. For both equations, we will present wellposedness theorems 
to the higher regularity. 

2.3.1. A-Stokes Equation. Let us consider the stationary Navier-Stokes problem. 

• Sa{p,u) = F in fl 

\/a • u = G in fl 

Sa{p-iU)M = H on S 

u = on S(, 

Since this problem is stationary, we will temporarily ignore the time dependence of rj^A, etc. 
Before discussing the higher regularity result for this equation, we should first define the weak formulation. 
Our method is quite standard; we will introduce p after first solving a pressureless problem. Suppose 
F G VF*, G G and H e H-'^/'^{T,). We say {u,p) eWxH° \sa weak solution to 2.3.1 if • u = G 
a.e. in fi, and 

(2.3.2) ^(]D'^w,P^v)„o - (p, • w)„o = (F,w)„o - {H,v)ho{i:) 
for all V €W. 

Lemma 2.21. Suppose F G W* , G G and H G i?~^/^(I]), then there exists a unique weak solution 
{u,p) €W xH° to 2.3.1. 

Proof. By lemma 2.17, there exists a. u € W such that ■ u = F^. Naturally we can switch the 
unknowns to w = u — u such that in the weak formulation w is such that Va ■ w — and satisfies 

(2.3.3) ^{BAW,BAv)no ~{p,Va- v)^^ = -^(D^u,D^t;)^o + {F,v)u'^ - {H,v)ho(i:) 
for all V eW. 

To solve this problem, we may restrict our test function to w G A such that the pressure term vanishes. 
A direct application of Riesz representation theorem to the Hilbert space whose inner product is defined 
as (u, v) = (lO^u, IiAv)-H° provides a unique w & X such that 

(2.3.4) i(P^u;,P^u)^o = -i(P^u,P^?;)„o + {F,v)ho - {H,v)no(j:) 
for all V G X. 

In order to introduce the pressure, we can define A G W* as the difference of the left and right hand sides 
in (2.3.4). So A(t;) = for all v G X. Then by proposition 2.20, there exists a unique p G H° satisfying 
{p, Va ■ v)u° = '^(^) for all v G W, which is equivalent to (2.3.3). □ 

The regularity gain available for solution to (2.3.1) is limited by the regularity of the coefficients of 
the operator A^, and V^-, and hence by the regularity of r]. In the next lemma, we will present 

some preliminary elliptic estimates. 

Lemma 2.22. Suppose that rj G ff''+^/^(E) for k > 3 such that the map $ defined in (1.3.3) is a G^ 
diffeomorphism of n to il' = $(11). If F e H°{n), G G H'^{^) and H G ffV2(s), then the equation 
(2.3.1) admits a unique strong solution [u.p) G X , i.e. {u,p) satisfies (2.3.1) in strong sense. 
Moreover, for r = 2, . . . , fc — 1 we have the estimate 

(2.3.5) \\u\\h. + |b||^._. < G{r,) \\F\\^^.. + ||G||^._, + ||^._3/.(j,) ) 

whenever the right hand side is finite, where G{ri) is a constant depending on ||??||iifi=+i/2(£;). 

Proof. This lemma is exactly the same as lemma 3.6 in [1], so we omit the proof here. □ 
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Notice that the estimate (2.3.5) can only go up to /c — 1 order, which does not fnlly satisfy onr 
requirement. Hence, in the following we will employ approximating argument to improve this estimate. 
For clarity, we divide it into two steps. In the next lemma, we first prove that the constant can actually 
only depend on the initial free surface. 

Lemma 2.23. Let k > 6 be an integer and suppose that rj G 7?''+^/^(E) and rjo e Then 
there exists eo > such that if ||?7 — 'no\\Hk-3/2^Y,^ < eo, solution to 2.3.1 satisfies 



(2.3.6) 



II^IIh^ + Mm-^ < C(r?o)( lli^lln^-. + + \\H\\ 



for r = 2,...,fc — I, whenever the right hand side is finite, where C{rio) is a constant depending on 

II'?o|Ih»»+i/2(E)- 

Proof. Based on lemma 2.22, we have the estimate 

(2.3.7) \\u\\h. + < C{v) \\F\\h.^. + ||G||^._, + ||i/||^._3/.(5:) ) 

for r = 2, . . . , A: — 1, whenever the right hand side is finite, where C{r]) is a constant depending on 

Define ^ = r? - ??o, then ^ e /f'=+V2(S). 
Let us denote and Ao of quantities in terms of r]o. We rewrite the equation 2.3.1 as a perturbation 
of initial status 

VA„-SAoiP,u) = F + F° m n 

Vao • m = G + go in n 

Sa„{p,u)J^o = H + H° on E 

u = on S(, 



(2.3.^ 



where 

F° = "^Ao-A ■ Sa{P, u) + Wao ■ Sao-aIP, u) 

(2.3.9) G" = V^„_^ • u 

i/o = Sa„ (p, u)iMo - AT) + SAo-AiP, u) 

Suppose that ||^||h'=-3/2(e) < 1, which implies ||^||^fe-3/2(E) < U\\hI'-^/^{t,) < 1 for any I > 1. A 
straightforward calculation reveals that 

||i^°||j;^.-2 < C(l + ll%llHfe+l/2(E))^ ||^||h^-3/2(s) {\\u\\h- + IbllHr-l) 

(2.3.10) < G(l + ||7?o|Ih^+i/2(j:))' UWh^-^/hi:) (II^IIhO 

||-H'||Hr-3/2(s) < C{1 + ||?70||Hfc + l/2(s))^ ll^llHfe-3/2(E) {11^^ + IblliJ'-l) 

for r = 2, . . . , fc — 1. 

Since the initial surface function % satisfies all the requirement of lemma 2.22, we arrive at the estimate 
that for r = 2, . . . , — 1 

(2.3.11) \\u\\„. + |b||^._, < C{r,o) ( \\F + F°\\^^_, + \\G + G°\\^^_, + \\H + H°\\^^_,,,^^^ ) 

where G(ryo) is a constant depending on ||r?o||jji=+i/2(s)- Combining all above, we will have 
(2.3.12) 

IMm + Ibll/f-i ^ 

C{V0) \\F\\Hr-2 + IIGIInr-l + ||/f||jj.-3/2(s) ^ + C(r?o)(l + ||??0||Hfc+l/2(E))* IICIIh^-3/2(S) {M Rr + Ibllffr-l) 

So if 

(2.3.13) U\\h--3/2(^) < min<{ ^ 



,._3/2(^) ^ \ 2 ' 4GG(,o)(l+b„||^...,.(,))^ 
we can absorb the extra term in right hand side into left hand side and get a succinct form 



(2.3.14) \\u\\j,r + Mnr-. < C{vo) ||F||h^-2 + ||G||^._, + ||//||h.-3/2(j:) 

for r = 2, . . . , fc - 1. □ 
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Note that the above lemma only concerns about regularity up to k — I and wc actually need two more 
order. Then the next result allows us to achieve this with a bootstrapping argument. 

Proposition 2.24. Let k > 6 be an integer, and suppose that ij G as well as rjo G 

satisfying \\r) — r?o||/f<!+i/2(x;) ^ eo- Then solution to 2.3.1 satisfies 

(2.3.15) \\U\\JJ. + WpW^r-. < C(%) ( \\F\\Hr-2 + ||G||^._. + ||iJ||H.-3/2(E) ) 

for r = 2, . . . ,k + 1, whenever the right hand side is finite, where C{r]o) is a constant depending on 

ll'?o||_H-f=+i/2(5]). 

Proof. If r < fc — 1, then this is just the conclusion of lemma 2.23, so our main aim is to gain two 
more regularity here for r = k and r = k + 1. In the following, we first define an approximate sequence 

for rj. In the case that S = R^, wc let p E C^{R^) be such that supp{p) C B(0,2) and p = 1 for 
B{0. 1). For m E N, define 77™ by r)'"(^) = p(£,/m)i)(S^) where ' denotes the Fourier transform. For each 
m, 77'" G H^{T,) for arbitrary j > and also 77™ ^ 77 in i?'^+^''^(S) as 7n — > 00. In the periodic case, 
we define 17"* by throwing away the higher frequencies: 7}™(7i) = for \n\ > m. Then 77™ has the same 
convergence property as above. Let A"^ and A/"™ be defined in terms of 77™. 

Consider the problem (2.3.1) with A and Af replaced by A"^ and A/"™. Since rj"^ G i/*=+5/2^ can apply 
lemma (2.22) to deduce the existence of (u"*,p'") that solves 



(2.3.16) 



and such that 



v^™ • Sa-- (p™, u") = f in n 

Va-u"^ = G in fl 

SA'^ip"\u"')Afo = H on S 

= on Sb 



for r = 2, . . . , k + 1. We can rewrite above equation in the following shape, as long as we split the EP^-nu" 
term. 



(2.3.18) 



-A^^w™ + V^mp™ = F + V^mG in fl 

V^™ • = G in n 

(p™7 - ©^mu'")A^™ = H on Y, 

= on Sb 



In the following, we will prove an improved estimate for (u'",p'") in terms of ||j?"*||ijfe+i/2(j])- We divide 
the proof into several steps. 

Step 1: Preliminaries 

To abuse the notation, within this bootstrapping procedure, we always use {u,p, 77) instead of {u"^,p"^, 77™) 
to make the expression succinct, but in fact they should be understood as the approximate sequence. 
Also it is easy to see the term V^mG will not affect the shape of estimate because |lV^mG||^fe_i < 
||77™|| j:^*;+i/2(x;) Hence, we still write F here to indicate the forcing term in first equation. 

We write explicitly each terms in above equations, which will be hired in the following. 
(2.3.19) 

diiui + d22Ui + {1 + A^ + B^)K^d33Ui - 2AKdi3Ui - 2BKd23Ui 

+{AKd3{AK) + BKdaiBK) - di{AK) - d2{BK) + Kd3K)dsUi + dip - AKdsp = Fi 

(2.3.20) 

diiU2 + d22U2 + (1 + ^2 + B^)K^d33U2 - 2AKdi3U2 - 2BKd23U2 

+ {AKd3{AK) + BKd3{BK) - di(AK) - d2{BK) + Kd3K)d3U2 + d2P - BKdsp = F2 

O'iOU + d22U3 + (1 + + B^)K^d33U3 ~~ 2AKdi3U3 " 2BKd23U3 

^ ^ ^ ^ +{AKd3{AK) + BKdaiBK) - di{AK) - d2{BK) + Kd3K)d3U3 + Kd3P = F3 

(2.3.22) diui - AKdaui + d2U2 - BKd3U2 + Kdsus = G 

where for all above. A, B and K should be understood in terms of rj"^. For convenience, we define 

(2.3.23) Z = G(r7o)P(r?) ( + ||G||^, + ||/f||^.-v2(E) ) 
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(2.3.25) 



where C{r]o) is a constant depending on ||?7o||jiffc+i/2(£;) and P{ri) is the polynomial of ||??||jiffc+i/2(£;). 
Step 2: r = k case 

For k — 1 order elliptic estimate, we have 

(2.3.24) Ml,., + Ml,., < Civo) ( ml,-s + ||G||h.-2 + ||/f||^._a/2(s) ) < ^ 

By lemma 2.23, the bounding constant C{r]o) only depend on ||??o||jjfc+i/2 • 
For i = 1,2, since (9,u, dip) satisfies the equation 

-AA{diu)+YA{diP) = F in fl 

Wa ■ [diu) = G in Q, 

{{dip)I - nA{diu))N = H on T. 

diU = on Sf, 

where 

F = diF + \7o,A ■ V^u + • \7o,AU - Vo.aP 
G = diG- Va^A ■ u 
H = diH- {pi - I}Au)diAf + Be,AuJ^ 
Employing fc — 1 order elliptic estimate, we have 

(2.3.26) \\diu\\l,., + \\dip\\l,., < C{vo) ( \\F\\1,-s + \\G\\1,., + ||^||?f.-a/2(s) ) 

Since except for the derivatives of F, G and H, all the other terms on the right hand side has the form 
\\A ■ B\\^r, 'm which A = d°'fj and B = d^u or d^p. These kinds of estimates can be achieved by lemma 
A.l. Because this lemma will be repeated used in the following estimates, we will not mention it every 
time and all of the following estimate can be derived in the same fashion. Hence, we have the forcing 
estimate 

||-^||jJfe-3 "I" ll^lljffc-2 "I" ||-^|| JJfe-5/2 

(2.3.27) < + + ||F||^.-3/2(^) + C{n,)P{n) ( + |b||^._2 ) 

<z 

In detail, this means 

(2.3.28) \\diu\\l,., + ||5ip||^._2 + \\d2u\\lu., + ||52p||^.-2 < Z 

Hence, most parts in + ||p|| jjie-i has been covered by this estimate, except those with highest order 

derivative of d^, ■ 

Multiplying A to (2.3.21) and adding it to (2.3.19) will eliminate the d^p term and get 

{diiui + Adiiu^i) + {d22Ui + Ad22U3) + (1 + + B^)K^{d33Ui + Adaau^) - 2AK{disUt + Ad^^u^) 
-2BK(d2?.u, + Ad2sUi) + {AKdi{AK) + BKd^iBK) - d^{AK) - d2{BK) + KdiK){d^u^ + Adim) 

+dip = Fi+ AF3 

Then taking derivative d^~^ on both sides and focus in the term (1 + + B^)K^[d^ui + Ad^us), the 
estimate of all the other terms in norm implies that 

(2.3.29) \\d^ui + Adlu3\\]^„ < Z 
Similarly, we have 

(2.3.30) \\d^U2 + Bd^usWljo < 2 
Rearrange the terms in (2.3.22), we get 

(2.3.31) K{1 + A'^+ B'^)d3U3 = G- dim - d2U2 + AK{d3Ui + Adsm) + BK{d3U2 + Bd3U3) 

Taking derivative ^3"^ on both sides, focusing in the term K{1 + A? + 5^)93^3 employing all of the 
estimate we have known, we can show 

(2.3.32) ||9|«3||Lo<^ 
Combining (2.3.29), (2.3.30) and (2.3.32), it is easy to see we can get 

(2.3.33) \\d^uif^, + \\dlu2f^,<Z 
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Plugging this to (2.3.21) and taking derivative ^3 ^ on both sides, we get 

|2 



1^3-^1 



(2.3.34) 

Combining this with all above estimate, we have proved 
(2.3.35) + 
Therefore, we have proved the case r = k. 

Step 3: r = k + 1 case: first loop 

For i,j = 1,2, since {dijU,dijp) satisfies the equation 



< :7 



-1 <^ 



(2.3.36) 



-A^(ay?i) + V_A{dijp) = F in fl 

■ (dij u) = G in fl 

{{d,jp)I-3A{d^ju)W = H on E 

(diju) = on S(, 



where 

F = dijF + Voi.A ■ ^Au + Va ■ "^d^.AU + VdiA ■ "^djAU + VdjA ■ Va.^u + VdiA ■ V^(5ju) 

+Va^^ • VA{diu) + Va ■ Vo^A{dju) + Va ■ ^d,A(.diu) - Vd^AP - ^diA{djp) - Vd^AidiP) 

G = dijG -Vg.^A ■ U-^diA ■ {dju) -S/g-A ■ (diu) 

H = d^H - {pi - nAu)d^jN - {{djp)I - BAidju) - Bo^Au)diM - {{dip)I - D^(5i«) - I])o,Au)djM 

+ {^d,,AU + Bg^AidjU) + na^A{^^u))N 

Employing k — 1 order elliptic estimate, we have 

\\diM\H>^-^ + \\di3P\\]i>^-^ P + G + H <Z 



(2.3.37) 



where the forcing estimate can be taken in a similar argument as in r = A; case. 

In detail, this is actually 



(2.3.38) \\diiu\ 



Hk-i 



\\di2u\\ 



Hk-1 



\dllp\\Hk 



\\di2P\\ 



Hk-2 



\\d22p\\Hk 



< Z 



Step 4: r = k + 1 case: second loop 

Similar to r = fc case argument, for i = 1, 2, taking derivative d^^^di on both sides of (2.3.29) and focus 
in the term d^diUi + Ad^dtu^, we get 

(2.3.39) \\dldiUi + AdldiUsW^^o < Z 
Similarly, we have 

(2.3.40) \dldiU2 + B5|9,«3||^o < Z 

Plugging in this result to (2.3.22) and taking derivative d\^^di on both sides, we will get 

(2.3.41) \d\d,u^f^,<Z 
An easy estimate for these three terms implies 

(2.3.42) \dld,uf^,<Z 

Plugging this to (2.3.21) and taking derivative d^^^di on both sides, we get 

(2.3.43) \\dt^d,pf^, < Z 

Combining all above estimate with estimate of previous steps, wc have shown 

(2.3.44) \\^3^ufH,-^ + ||53ip||^.-2 < Z 

In detail, this is actually 

(2.3.45) \\di3u\\l,., + \\d23u\\l,., + ||ai3p||^._. + \\d23p\\l,-. ^ Z 
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Step 5: r = k + I case: third loop 

Again we use the same trick as above. Taking derivative d^~^ on both sides of (2.3.29) and focusing in 
the term Sg+^ui + A^g+^ua, we can bound Sg+^ui - Ad^^^us. Similarly, we control 53+^^2 - Bd'^^^u^.. 
Plugging in this result to (2.3.22) and taking derivative ^3"^ on both sides, we can estimate d'^'^^u^. 
Then we have 

(2.3.46) \\dl+'u\^^,<Z 

Plugging this to (2.3.21) and taking derivative ^3"^ on both sides, we get 

(2.3.47) \\dtp\^^,<Z 
Combining this will all above estimate, we get 

(2.3.48) 11^33^x11^.-1 + ||533p||^.-2 < Z 



Step 6: r = fc + 1 case: conclusion 

To synthesize, (2.3.38), (2.3.45) and (2.3.48) imply that all the second order derivative of m in fc — 1 norm 
and p in fc — 2 norm is controlled, so we naturally have the estimate 

(2.3.49) ll^^llUi + Ibll^^ ^ C{i^^)P{v) ( II J^ll^.-i + l|G||^. + ||if ||^.-V2 

This is just what wc desired. 

Now let us go back to original notation, since P{-) is a fixed polynomial, this gives an estimate 



(2.3.50) 



< C(r7o)P(r/) \\F\\\,^, + ||G||^. + \\H\\\,^,,.^^^ ^ 

< C(r7o)P(r?o)( ||f ll^.-i + ||G||^. + ||if ||^.-v2(s) 
<C{vo)(\\F\\l.-r + \\G\\l, + \\H\\ 



where C(?7o) depends on ||??o||//fc+i/2(x;)- 

This bound implies that the sequence (u™,p™) is uniformly bounded in x H'^, so we can extract 

weakly convergent subsequence u"^ u'^ and p™ . 

In the second equation of (2.3.16), we multiplied both sides by J"^w for w S to see that 

/ GwJ™ = / (V^m • u'")wJ™ = - / • (V^mw) J™ [ vP . {Vaw)J = [ {Va- u°)wJ 
Jn Jn Jn Jn Jn 

(2.3.51) 

which implies Va ■ vP = G. Then we multiply the first equation by w J™ for w e and integrate by 
parts to see that 

(2.3.52) i / P^n-u™ : D^^wJ™ - / p^V^m • wJ™ = / F • w J" - f H ■ w 
2 Jo, Jn Jn Jn 

Passing to the limit m ^ 00, we deduce that 

(2.3.53) i / Bau° : BawJ - / p°Va ■ wJ = / F ■ wJ - I H-w 
2 Jn Jn Jn Js 

which reveals, upon integrating by parts again, (u°,p°) satisfies (2.3.1). Since {u,p) is the unique strong 
solution to (2.3.1), we have vP = u and p° = p. This, weakly lower semi-continuity and estimate (2.3.50) 
imply (2.3.54). □ 

Remark 2.25. The key part of this proposition is that as long as we can deduce 770 € H^+^/'^{T) and 
•q e i?*^+^/^(S), we can achieve u G 7?*^+^ and p G H^, which is the highest possible regularity we can 
expect. 
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Remark 2.26. By our notatum, rule, since C{riQ) depends on Qq and is given implicitly, we can take it 
as a universal constant, so the estimate may be rewritten as follows. 

(2.3.54) llt^ll^. + M^r-. < \\F\\^... + ||G||^._. + ||ff||H.-3/2(s) 

for r = 2, . . . , fc + 1 

2.3.2. A-Poisson Equation. We consider the elliptic problem 

{A^p = / in Q. 
p = g on S 

V^p ■ V — h on Sf, 

For the weak formulation, we suppose / € V* , g e 7J^/^(S) and h e i?~^/^(Sb). Let p G he 
an extension of g such that supp{p) C {— inf(&)/2 < xs < 0}. Then we can switch the unknowns to 
q = p — p. Hence, the weak formulation of (2.3.55) is 

(2.3.56) {VaQ, '^A<i>)n° = -{^aP, V^</')„o - (/, + {h, (P)_y2 

for all (j) V, where (•, ■)v denotes the dual pairing with V and (•, •)_i/2 denotes the dual pairing with 
i?~^/^(S(,). The existence and uniqueness of a solution to (2.3.56) is given by standard argument for 
elliptic equation. 

If / has a more specific fashion, we can rewrite the equation to accommodate the structure of /. Suppose 
the action of / on an element (p G V is given by 

(2.3.57) = {.fo,(^)no + (Fo,V^0)„o 

where (/o,-Fo) & H° x H° with ||/o||^o + Hi^oUffo = Then we rewrite (2.3.56) into 

(2.3.58) {VaP + Fo, Va<P)h° = "(/o, + {h, </.)_i/2 
Hence, it is possible to say p is a weak solution to equation 

f • (V^p + Fo) = /o 

(2.3.59) I P = g 

I (V^p + Fo)-i^ = /i 

This formulation will be used to construct the higher order initial conditions in later sections. 

Lemma 2.27. Suppose that rj € H'^~^^^'^{'E) for k > 3 such that the map ^ is a diffeomorphism of 
Q to 0,' = $(0). If f G H^,g e H^/^ and h € H^^'^, then the equation 2.3.55 admits a unique strong 
solution p e H'^. Moreover, for r = 2, . . . , — 1, we have the estimate 

(2.3.60) IIpII^. < C{ri) ||/||^._. + WgWnr-u^,:) + Mnr-v^,:,) ) 
whenever the right hand side is finite, where C{ri) is a constant depending on \\v\\h''+^/^{^)- 

Proof. This is exactly the same as lemma 3.8 in [1], so we omit the proof here. □ 

Next, we will prove the bounding constant for the estimate can actually only depend on the initial 
surface. Since we do not need optimal regularity result for this equation, we do not need the bootstrapping 
argument now. 

Proposition 2.28. Let k > 6 he an in teger suppose that t] e H''+^/'^{E) and % e iJ'=+i/2(S). There 
exists Co > such that if \\rj — ryoU < eo, then solution to 2.3.55 satisfies 

(2.3.61) lblli/.<C(%)(||/||H^-. + NI 

for r = 2, . . . ,k — 1, whenever the right hand side is finite, where C{r]o) is a constant depending on 

ll'?o|lHfc + i/2(H) • 

Proof. The proof is similar to lemma 2.23. We rewrite the problem as a perturbation of the initial status. 

r A^„p™ = / + /'» in n 

(2.3.62) } p^ = g + g"' on 1: 

[ VaoP"" ■v = h + h"' on Eb 

The constant in this elliptic estimate only depends on the initial free surface. We may estimate f™, g™ 
and /i™ in terms of p"*, rjo and ^"^ = 77™ — ryo- The smallness of ^"^ implies that we can absorb these 
terms into left hand side, which is actually (2.3.63). □ 
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Remark 2.29. Similarly, we can also take this constant C(ryo) as a universal constant and rewrite the 
estimate as follows. 

(2.3.63) Ibll^. < ||/||^._. + + ll/j|lHr-3/.(s,) 

for r = 2, . . . , fc — 1. 

2.4. Linear Navier-Stokes Estimates. 

dfU — Aj[u + Vytu = F in fl 
■ u — in n 



on S 

u = on 'Eb 

In this section, we will study the linear Navier-Stokes problem (2.4.1). Following the path of I. Tice and 
Y. Guo in [1] . Wc will employ two notions of solution: weak and strong. Then we prove the wcUposcdncss 
theorem from lower regularity to higher regularity. To note that, most parts of the results and proofs 
here are identical to section 4 of [1] , so we won's give all the details. The only difference is the constants 
related to free surface since our surface has less restrictions. 

2.4.1. Weak Solution and Strong Solution. The weak formulation of linear problem (2.4.1) is as follows. 

(2.4.2) {dtu,v)c2-H° + ^{^AU,'^Av)c''H° - {P,'^A ■ v)c'^n'' = {F,v)c''n° - {H,v)c2nocE) 

Definition 2.30. Suppose that uq G H", F e W* and H e L^{[0,T]; H-^/'^{T,)). If there exists a pair 
{u,p) achieving the initial data uq and satisfies u € W, p G L^{[0,T]; and dtU G W* , such that the 
(2.4-2) holds for any v <E W, we call it a weak solution. 

If further restricting the test function v G X, we have a pressureless weak formulation. 

(2.4.3) {dtU, v)c^u° + ^{^AU, ID)^u)£2^o = {F, v)^^-^^ - {H, w)£2^o(s) 

Definition 2.31. Suppose that uo e y{0), F e X* and H G L'^{[0,T]; H-'^/'^CE)). If there exists a 
function u achieving the initial data uq and satisfies u G X and dtU € X* , such that the (2.4.3) holds 

for any v Cz X , we call it a pressureless weak solution. 

Remark 2.32. It is noticeable that in fi, W* C X* and \\u\\j^, < || 

^llw*' ^ ^ ' have 

Since our main aim is to prove higher regularity of the linear problem, so the weak solution is a natural 
byproduct of the proof for strong solution in next subsection. Hence, we will not give the existence proof 

here and only present the uniqueness. 

Lemma 2.33. Suppose u is a pressureless weak solution of (2.4.1). Then for a.e. t € [0,T], 
(2.4.4) 



Also 



^t; I J\u{t)\' + - I J\B^u\' = - m\uo\' + - \u\'d,J+ J{F-u)~ I jH-u 

'a JQ ^ Jo JQ Jo Jn 



(2.4.5) ll^^lli^^o + Ml,^, < Coiv) ( hollno + lli^ll^- + ml^H-^/Hi:) ) 
where 

(2.4.6) Coiv) = P(||»?|IhV2(e) + IISt^Jllnv^fE)) exp (TP{\\dtv\\H^;,^^))) 
Hence, if a pressureless weak solution exists, then it is unique. 

Proof. The proof is almost the same as that of lemma 4.1 and proposition 4.2 in [1]. The only difference 
is that we should replace the lemma 2.1 in [1] used in that proof by our lemma 2.11 now. □ 

Next, we define the strong solution. 

Definition 2.34. Suppose that 

uoGH'^n X{0) 

(2.4.7) F G L'^{[0,T];H^)nC°{[Q,T];H°) dtF £ X* 

H e L2([o,T];if3/2(5])) nC0([0,T];ifV2(s)) dtH e ^^([O, T]; /f-V2(S)) 
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// there exists a pair {u, p) achieving the initial data uq and satisfies 

ueLH[0,T];H^)nC^{[0,T];H^)nA: d,u e L^{[0,T]; H') n C"{[0,T]; H'^) 
^ ^ dfueX* pe L^{[i),T]:H^)nC%[0,T];H^) 

such that they satisfies (2.4-1) in the strong sense, we call it a strong solution. 

2.4.2. Lower Regularity Theorem. 

Theorem 2.35. Assume the initial data and forcing terms satisfies the condition (2.4-7). Define the 
divergence-A preserving operator DfU by 

(2.4.9) Dtu = dtu-Ru for R = dtMM-^ 

where M is defined in (2.2.23). Furthermore, define the orthogonal projection onto the tangent space of 

the surface {x^ = ??o} according to 

(2.4.10) Iio{v) ^ V - {v No)MuWol^ for A^o = (-^1770, -92%, 1) 
Suppose the initial data satisfy the compatible condition 

(2.4.11) no(ir(0) + B^oUoA/'o) = 

Then the equation (2-4-1) admits a unique strong solution {u,p) which satisfies the estimate 
(2.4.12) 

\H\I2H3 + \\dtu\\l2Hi + \\dfu\\^, + \\u\\l^^2 + \\dtu\\l^ijo + \\p\\l2h^ + \\p\\l<^m 
where 

(2.4.13) L(?7) = P{K{v)) exp (TP{K{r])) 
for 

(2.4.14) if (r?) = ^sup^ h(t)||^o/.(s) + \\dtrim 

The initial pressure p{Q) € H^, is determined in terms of uq, rjo, F{0) and H{0) as the weak solution to 
( Vx„ • {Va„p{0) - F{0)) = -VA„iRiO)uo) in n 

(2.4.15) <^ piO) = {H{0)+BAoUof^o)-J^o\^ol^ on E 
[ {WaoP{0)-F{Q))-^ = ^AoUo-t^ on Sb 

in the sense of (2.3.59). 

Also, Dtu{Q) = dtu{0) — R{0)uo satisfies 

(2.4.16) DM0) = AaoUo - VaoP{0) + F{0) - R{0)uq G 3^(0). 
Moreover, {Dtu, dtp) satisfies 

( dt{Dtu)-l^A{Dtu) + VA{dtp) = DtF + G^ in Q 

/,4 17X J ^A-{Dtv)=0 in ft 

^ ' I SA{dtp,Dtu)M = dtH + G" on S 

[ Dtu = on Tif, 

in the weak sense of (2.4-3), where 

= -{R + dtJK)I^AU-dtRu + {dtJK + R + R^)VAP + '^A-{^A{Ru)-mAU + ^a,Au) 
G" = Ba {Ru)J\f - {pi - BAu)dtJ\f + BdtA-uJ^ 

Proof- Since this theorem is almost identical to theorem 4.3 in [1], we will only point out the differences 
here without giving details. In order for contrast, we use the same notation and statement as [1] here. 

(1) The only difference in the assumption part is that in theorem 4.3 of [1], K{r]) is sufficiently small, 
but our K{ri) can be arbitrary. Hence, in [1], we can always bound P{K{r])) < 1 + K{'n)^ but 
now we have to keep this term as a polynomial (this is the only difference of L{ri) between these 
two theorems). 

(2) It is noticeable that each lemma used in proving theorem 4.3 of [1] has been recovered or slightly 

modified in our preliminary section. To note that these kind of modification merely involves the 
bounding constant due to the arbitrary K{r]), but does not change the conclusion and the shape 
of estimate at all. Therefore, all these modification will only contribute to P{K{r])) term. 
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(3) In elliptic estimate, we now use proposition 2.24 to replace proposition 3.7 in [1], which gives 

exactly the same estimates. 

Based on all above, we can easily repeat the proving process in [1] with no more work, so our result 
naturally follows. □ 

2.4.3. Initial Data and Compatible Condition. We first define some quantities related to the initial data: 



(2.4.18) 
(2.4.19) 

(2.4.20) 
(2.4.21) 

(2.4.22) 



Ho{u,p) 
/Co(uo) 

^o(r?) 

no{F,H) 



\\uo\\h2N 



J=0 



fj2N-2j-l 



N 



||r?(0)||^.™(^)+E|^''?(0) 

i=i 

||??0||//2N + l/2(2) 



^2N-23+3/2(J]) 



Af-1 
j=0 



N-1 



^2iV-2j-2 



E FtH{o) 



fJ2N-2j-3/2(Y,) 



Furthermore, we need to define mappings Qi on O and Q2 on E. 
(2.4.23) 

g^{v, q) = -{R + dtJK)AAV - dtRv + {dtJK + R + R^)\l a(1 + • {^a{Rv) - R'^av + ^a,Av) 

g'^{v,q) = Oa{Rv)N - {ql - OAv)dtN + Oo^AvM 
The mappings above allow us to define the arbitrary order forcing terms. For j = 1, . . . , N , we have 
F° = F 

(2.4.24) ^ ^ ^pip) ^ ^ ^.-^1 j^i^g^ (Dr'-'u, dr'-'p) 

w = DtW-^ + g\Dr\, dr'p) = DiH + ECo Dig2{Di-'-\, df-'-'p) 

Finally, we define the general quantities we need to estimate as follows. 



N 



(2.4.25) £{u,p) 

(2.4.26) V{u,p) 

(2.4.27) IC{u,p) 



2 



N-1 



2. + E i^*^'^ 



j=0 



N 



j;,OOJJ2iV-2j-l 

N-1 



3=0 



+ \K+'u\\" + y Pip 

£{u,p)+V{u,p) 



I^2IJ2N-2j 



N 



(2.4.28) 

(2.4.29) 
(2.4.30) 



^(»?) = ll»7llLooH2iv+i/2(s) + E 



P(r?) 

/C(r,) 



L«>/f2N-2j + 3/2(2) 
N+l 



II^IIl2H-2N + 1/2(s) + ||9tr/||^2H-2N-l/2(s) + E 

f(r?)+2?(77) 



i2JJ2JV-2j+5/2(2) 



(2.4.31) 1C{F,H) 



N-l 

E l^t^ 
3=0 

N-l 

+ E 1^*^ 
3=0 



2 

L2H^N- 



N 



N-l 

+ y 

2j-2 ^-^ II ^ 



3=0 



j=0 
2 

j;,OOJj2]V-23-3/2(J]) 



J,2JJ2N-2j-l/2(s) 
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Before discussing the higher regularity, wc first define the compatible condition of initial data. For 
j = 0, . . . , Ai" — 1, we say that the j*^ compatible condition is satisfied if 



^ ' ' ' I no(//J(o)+D^„D>( 



Uo{H^{0)+BA,Diu{0)^^o) = 

Then we may construct the initial condition of higher order satisfying the compatible condition. Since 

this section is almost identical to section 4.3 in Y. Guo and I. Tice's paper (see pp 34-35 of [1]), we omit 
the details here. The core idea is to construct the initial data iteratively from lower to higher order. 
Then we have the estimate. 

(2.4.33) no{u,p) < P{Ho{v)) ( IkoHn^Jv + noiF, H)^ 
2.4.4. Higher Regularity Theorem. 

Theorem 2.36. Suppose the initial data uq G H^^ , rjQ € ff'^^+"'^/^(S) and the initial condition is 
constructed as above to satisfy the j*'' compatible condition for j = 0,1,...,A^ — 1, with the forcing 
function K,{F,H) + ICo{F, H) < oo. Then there exists a universal constant To{ri) > depending on C{t]) 
as defined in the following, such that if < T < To{r]), then there exists a unique strong solution {u,p) 
to the equation (2.4-1) on [0,T] such that 

5/wGL2([o,T];i/2JV-2i+i)p^O([o,T];i/2^-2j) j = 0,...,N 

(2.4.34) a/pG L2([0,T];if2iv-2j-)pC'0([0,r];//2iv-2J-i) for j = 0,...,N-l 

Also the pair {u,p) satisfies the estimate 

(2.4.35) lC{u,p) < £(77) ||wo|Ih2« + no{F, H) + /C(F, H) 
where 

(2.4.36) C{r,) = P{K{ri)) exp f TP(/C(»7)) 



(2.4.37) 



Furthermore, the pair {Dlu,dfp) satisfies the equation 

dt{Diu)-Aj^{Diu) + VA{dip) = F^ in Q 

Va ■ {D{u) =0 in n 

SA{dip, Dlu)N = Hi on S 

Din = on Sfe 

in the strong sense for j = 0, . . . ,N — 1 and in the weak sense for j = N. 

Proof. Similar to lower regularity case, this theorem is almost identical to theorem 4.8 in [1]. Due to 

the same reason that our free surface 77 can be arbitrary, the only difference here is that we keep the 
polynomial of /C(ry) and T-Lo^r]) in the final estimate which cannot be further simplified. Hence, our result 
easily follows. □ 

2.5. Transport Equation. In this section, we will prove the wellposedness of the transport equation 

dtri + uidiTj + ^2^2?? = U3 on E 
rj{t = 0) = ?7o on S 

Throughout this section, we can always assume u is known and satisfies some bounded properties, which 
will be specified in the following. 

2.5.1. Wellposedness of Transport Equation. Define the quantity 



(2.5.1) 



(2.5.2) Q(u) = y^\\diu ' + V \\dl 
Obviously, we have the relation Q{u) < lC{u,p). 

Theorem 2.37. Suppose that ICo(jio) < 00 and Q{u) < 00. Then the problem (2.5.1) admits a unique 
solution T] satisfying IC{r]) < 00 and achieving the initial data dfr]{0) for j = 0, . . . ,N . Moreover, there 
exists a < T{u) < 1, depending on Q{u), such that ifO<T<T, then we have the estimate 

(2.5.3) /C(r?) < P(/Co(%)) + P(S(«)) 
where P(-) is a polynomial satisfying P{Q) = 0. 
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Proof. We divide the proof into several steps. 
Step 1: Solvability of the equation: 

Similar to proof of theorem 5.4 in [1], based on proposition 2.1 in [5], since u e L^([0, T]; 

then the equation (2.5.1) admits a unique solution rj S L°° ([0,T]; H^^~^^/^(T,)) achieving the initial data 

v{0) = m- 

Step 2: Estimate of f (?y): 
By lemma B.l, we have 

(2-5.4) ||r?||^„,jj2iv+i/2(s) <exp(^C^ ||u(t)||^2Ar+i/2(s) dt^ (^V^Vo) + ||M3(i)|lH2iv+i/2(s) c^* 
for C > 0. The Cauchy-Schwarz inequality implies 

(2.5.5) C \Ht)\\H2M+u2(^E)dt<C \\u{t)\\H2N+idt<Vf^/Q(JI) 
So if T< l/(2Q(w)), then 

(2.5.6) ^""^l*^/ ll«Wllff2«+i/2(s)rfi) <2 
Hence, we have 

(2.5.7) M\i </Co(r?o)+TQ(w) </Co(r?o) + l 
Based on the equation, we further have 

(2.5.8) \\dtr]\\j^ooH2N-l/2(^s) < ||u3|Il'=oH21V-1/2(s) + I|w|IlooJJ2JV-1/2(s) ll'?llz,c»Ji-2]V + l/2(J]) 

< ICo{Vo)Q{u) + Q{u)+TQ{uf 

< (iCo{vo) + i]q{u) 



The last inequality is based on the choice of T as above. The solution rj is temporally differentiable to 
get the equation 

(2.5.9) dt{dtri) + uidi{dtr]) + U2d2{dtri) = dtu^ - dtuidirj - dtU2d2T] 
So wc have the estimate 

(2.5.10) ll^t^??!! j;,oojy2N-5/2(2) ~ 1 1 " 1 1 L°° H^N - 5/2 (J]) 1 1 9t 1 1 ^oo ^2iV - 3/2 (S) + 1 1 " 1 1 L=° "6/2 (S) 

+ ll^t^llL°°i/2]V-5/2(s) ||?7ll L^_ff2W-3/2(2) 



In a similar fashion, we can achieve the estimate for dfrj as j = 1,. . . ,N. 

(2.5.11) dfj] ' < f/Co(r?o) + Q{u) 

^ ' t ' j;,ooJj2]V-23+3/2(J]) ~ ^ "V'"/ / 

To sum up, we have 

(2.5.12) S{rj)<{l+ICo{vo))P{Q{u)) 



Step 3: Estimate of V{r]): 

Similar to above argument, by differentiating on both sides and utilize the transport equation iteratively, 

we can show 

(2.5.13) 'D{v)<{l + ICoivo))PiQiu)) 
Naturally, summarizing all above, we have 

(2.5.14) /C(r?) = £{r,) + V{r,) < (1 + /Co(r?o))P(Q(«)) 

Then our result easily follows. □ 
Next, we introduce a lemma to describe the different between r] and r]o in a small time period. 
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Lemma 2.38. // Q{u) + /Co(??o) < oo, then for any c > 0, there exists a T'(e) > depending on Q{u) 
and /Co(r?o)j such that for any < T < T, we have the estimate 

(2.5.15) ||?7 - ??o|lioojy2jv+i/2(s) < e 

Proof. ^ = T] — rjo satisfies the transport equation 

(o K^R\ / + ■"i^iC + ■"2^2^ = U3- uidir]o - W2C^2% on E 

^^•^■^^^ \ m = 

By the transport estimate stated above, we have 



(2.5.17) 



m\ 



i;,ooj/2iV+l/2(s) 
T 



< exp C 



Ht)\\ 



^2JV + l/2(2) 



dt 



□ 



i|u3(^) - ui{t)dir]o - '«2(t)52?7o|l//2N+i/2(£) 

< TQ{u)IC„irjo) 
Hence, when T is small enough, our result naturally follows. 

2.5.2. Forcing Estimates. Now we need to estimate the forcing terms which appears on the right-hand 
side of the linear Navier-Stokes equation. The forcing terms is as follows. 

(2.5.18) F = dtfjbKdsu-u-VAU 

(2.5.19) H = rjAf 
Recall that we define the forcing quantities as follows. 

2 

j;^2//2W-2j-l 



(2.5.20) IC{F, H) = 



N-l 

3=0 



+ \\d^F 



|2 



N 
3=0 



j;^2//2W-2j-l/2(s) 



j=0 



(2.5.21) no{F,H) 



N-l 

El 

3=0 



diFiO) 



L°°//2JV-2j-2 



N-l 

+ E 

3=0 

N-l 

E" 

3=0 



\dlH 



dlH{Q) 



2,oo/^2N-2j-3/2(2) 



//2JV-2j-3/2(s) 



In the estimate of Navier-Stokes-transport system, we also need some other forcing quantities. 



(2.5.22) T{F, H) = 



(2.5.23) U{F, H) 



N-l 

= T.\ 

3=0 



WfF 



N-l 

E 1^*^ 

3=0 



^2jLf2JV-2j-l 



£,2//2N-2j-l 



+ Wf 



+ 



N 

El 

3=0 



L°°//2N-2j-l/2(s) 



Af-1 

E 1^*^ 

3=0 



i2//2JV-2j-l/2(x;) 



Theorem 2.39. The forcing terms satisfies the estimate 



(2.5.24) 
(2.5.25) 
(2.5.26) 



noiF,H) 
TiF,H) 



< 
< 
< 

< 



P{ic{v)) + P{Q{u)) 

P(/Co(??o))+P(/Co(wo)) 
PilCiv)) + PiQiu)) 

rfp(/C(r?))+P(Q(w)) 



(2.5.27) HiF,H) 
where P(-) is a polynomial with P{0) = 0. 

Proof. The estimates follow from simple but lengthy computation, involving standard argument, so we 
will only present a sketch of the proof here. 

After we expand all the terms with Leibniz rule and plug in the definition of F^ and , it is easy to see 
that the basic form of estimate is ||Xy||^2jyfc or HXyH^oo/fi. where X includes the terms only involving 
fj and Y includes the terms involving u, p, F or H. Employing lemma A.l, we can always bound 



l^lliocfffei \\y\\l2Hi'2, and also llXyH^ooH-e ^ 



II A"||^3o jffci ||y|lioo^fc2 • The principle of choice in these two forms for L H is not to exceed the order of 
derivative in the definition of right-hand sides in the estimates. A detailed estimate can show this goal 
can always be achieved. In (2.5.24), note the trivial bound 



(2.5.28) 



d^Fl 



< Wf 
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Also in (2.5.27), the key part is the appearance of bounding constant T. We first trivially bound it as 
follows. 

Af-l „ Af-1 



j=0 



2 

N-1 



2j-l /-^ II ^ 



i=o 

2 



Z,2H2W-2i-l/2(s) 
JV-1 



Finally, an application of inequality 2a6 < + will imply the required estimate. 



iooJf2N-2j-l/2(x;) 



□ 



Remark 2.40. The reason why we can get the constant T lies in that, in the nonlinear Navier-Stokes 
equation, the u in nonlinear terms actually has one less derivative than the linear part, which makes 
it available to use lemma A. 15. The appearance of T in (2.5.27) will play a key role in the following 

nonlinear iteration argument. 

2.6. Navier-Stokes-Transport System. In this section, we will study the nonlinear system 

dtu — dtfjbKd^u + u ■ js,u — A^u + V^p = in f2 

■ u — Q in f2 

Sa{p, u)Af = fjAf on S 

u = on Sb 
u{x, 0) = ^0(2;) 



on 



(2.6.1) 



diTj + uidiT] + U2d2r] — U3 

We will first show a revised version of construction of the initial data and then employ an iteration 

argument to show the wcUposedncss of this system. 

2.6.1. Initial Data and Compatible Condition. Before we turn into the nonlinear equation, we first need 
to determine the higher order initial condition. 

Define 

(2.6.2) /Co=/Co(Mo)+/Co(r?o) 
We say (uo,%) satisfies the A''*'' order compatible condition if 

r • (D^w(o)) = in n 

(2.6.3) I Diu{0) = on 

i Ilo{W{0)+-BAoDiu{0)N'o) = on S 

for j = 0, . . . , N — 1. Then wc can employ an iterative argument to construct the initial condition 
of higher order. We will not repeat the detailed process here (see section 5.2 of [1]) and only give an 
estimate for the constructed data. 

Theorem 2.41. Suppose {uq, r]o) satisfies /Co < 00. Let w(0), df r]{0) for j = 0, . . . ,N and dlp{Q) for 
j = 0, . . . ,N — 1 defined as we stated. Then 

(2.6.4) /Co < no{u,p) + Hoiv) < P{}Co) 
for a given polynomial P(-) with P{0) = 0. 

Proof. This is a natural corollary of our construction, so we omit the proof here. □ 

2.6.2. Construction of Iteration. For given initial data («o,%), by extension lemma A. 16, there exists 
w° defined in f2 x [0, 00) such that Q(m°) < /Co(tio) achieving the initial data to N*^ order. By solving 

transport equation with respect to theorem 2.37 implies that there exists 77" defined in f2 x [0,To) 
such that ICir]°) < (1 + /Co(%))P(Q(w")) < P{ICo) < 00. This is our start point. 

For any integer m > 1, define the approximate solution (m™,p™, rf ) on the existence interval [0, T^) by 
the following iteration. 

' 9tM™- A^^-iu^ + V^m-ijj™ = 5tr?'"-^6ii''"-^53U™-^ in n 
V^^-i • = in n 

(pmj _ D^™-iu™)j\^™-i = on S 



(2.6.5) 



= 



on 



on 
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where {u"^,ri"^) achieves the same initial data (uo-r^o). and A"^' . A/"'", A'"' are given in terms of t/™. 
This is only a formal definition of iteration. In the following theorems, we will finally prove that this 
approximate sequence can be defined for any m e N and the existence interval will not shrink to 
as m — >■ oo. 

2.6.3. Boundedness Theorem. 

Remark 2.42. Before we start to prove the boundedness result, it is useful to notice that, based on the 
linear estimate (2.4-35) and forcing estimate in lemma 2.39, this sequence can always he constructed and 
we can directly derive an estimate 

(2.6.6) K{vr+\p^+^) + lC{v"'+^) < PiQiu"") + /C(r?™) + /Co) 

when T is sufficiently small, where P{-) is a polynomial satisfying P{0) = 0. Since the initial data can 

be arbitrarily large, this estimate cannot meet our requirement. Hence, we have to go back to the energy 
structure and derive a stronger estimate. However, this result naturally implies a lemma, which will be 
used in the following: 
If 

(2.6.7) Q(u"-^) + /C(r/'"-i) < Z 
then 

(2.6.8) /C(u'",p'") + /C(r?™) < CP(/Co + Z) 

Theorem 2.43. Assume J° > 6 > and the initial data {uo,rio) satisfies the N*^ compatible condition. 
Then there exists a constant < Z < oo and < T < 1 depending on ICq, such that if < T < T and 
/Co < 00, then there exists an infinite sequence {u"^,p'^,ri"^)^^Q satisfying the iteration equation (2.6.5) 
within the existence interval [0, T) and the following properties. 

(1) The iteration sequence satisfies the estimate 

Q(u")+/C(rr) <Z 

for arbitrary m, where the temporal norm is taken with respect to T . 

(2) For any m, J'^it) > 6/2 forO<t<T. 

Proof. Let's denote the above two assertions related to m as statement P^. We use induction to prove 
this theorem. To note that in the following, we will extensively use the notation P(-), however, these 
polynomials should be understood as explicitly given, but not necessarily written out here. Also they 
can change from line to line. 

Step 1: Po case: 

This is only related to the initial data. Obviously, the construction of u° leads to Q(u°) < /Co. By 
transport estimate (2.37), we have K,{rf) < -P(/Co) for P(0) = 0. We can choose Z > P(/Co), so the first 
assertion is verified. 

Define = rf — rio the difference between free surface at later time and its initial data. Then the estimate 
in lemma 2.38 implies ||C°|Looh5/2 < TZ. Naturally, sup^gfo^y] | J0(i) - J°(0)| < ||^°||^^^5/2 < TZ. 
Thus if we take T < S/{2Z), then J'^(f) > 6/2. So the second assertion is also verified. In a similar 
fashion, we can verify the closeness assumption we made in linear Navier-Stokes equation, i.e. rj and r]o 
is close enough within [0, T] can also be verified. Hence. Po is true. 

In the following, we will assume Pm-i is true for m, > 1 and prove P„i is also true. As long as we can 
show this, by induction, P„ is valid for arbitrary n e N. Certainly, the induction hypothesis and above 
remark implies Q(u'"-i) + /C(7?'"-i) < Z and /C(u'",p'") + /C(r?'") < CP(/Co + Z). 

Step 2: Pm case: estimate of via energy structure 

By theorem 2.36, the pair {D^u"^,d^u™) satisfies the equation (2.4.37) in the weak sense, i.e. 



(2.6.9) 



' at(Dfw™)- A^n.-i(£)fO + V^m-i(5fp'") = in SI 
V^m-i • {D^tu"^) =0 in n 

(£»f w™, d^p"")}^ =H^ on 

,r 
t 



D^vJ^ = on Sfe 
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where and H'^ is given in terms of u™ and rj™ ^ . Hence, we have the standard weak formulation, 

i.e. for any test function tp G the fohowing holds 

(2.6.10) 



Therefore, when we plug in the test function i/; — D^u"^, we have the natural energy structure 
(2.6.11) i / j\D^u"^\' + l f I J|D^.-ii5f«-|' = 



J(0) iD^u'^iO) 



2 Jo 



N„,m I 



jpN . n^^m _ 



JO. 



Jn 



A preliminary estimate is as follows. 



N„,m\ 



LHS > \\D-tU 



\L°°H« 



RHS < P{ICo) + TZ\\D^u 



-^\/t7\\f^\\ II n^?/™!! -I- \/t II TT^II Wri^ii'^W 

-t-V-t^ll-f 11^2^0 ll-t^t 1* Ili^ffO + V J IP |lL°°ff-l/2(S) ll-^t " 11^2^1/2(2) 



< P{Ko) + TZ\\D^u 



+ VrllF- 



+ Vtz'^ \\d. 



+ y/T\\H 



li:,oojj-l/2(£;) 



+Vf\\D^vr\ 



for a polynomial P(0) = 0. Taking T < 1/(16^^) and absorbing the extra term on RHS into LHS implies 
(2.6.12) \\D^u^\\1^^, + \\D^u^\\1,^,<P{Ko) + Vt\\F''\\%^, + Vt\\H^\^^^^ 



Drop the tt 



N„,m I 



term and we derive further the estimate for d^u^ 



(2.6.13) ll^r^i^ll'.^i <P(/Co) + Vr||i^^| 



Vt\\h 



N\ 



+ II D^?;" - (9^7;™l 

Then we need to estimate each term on RHS. For the middle two terms, it suffices to show it is bounded, 
however, for the last term, we need a temporal constant T within the estimate, which can be done by 
lemma A. 15. Since these estimates is similar to the proof of lemma ??, we will not give the details here. 

■ N-l „ N-1 



\F 



N 



< P{Ko+Z)+T{F,H) 



2 

L2ii-2 



j=0 



2 

1,2 J/l 



-:f{f,h) 



lL~i?-i/2(E) 



(iv — 1 

< Pi]Co+Z)+J^{F,H) 



2 



N-1 

■ E P'tP" 

j=0 



2 



. N-1 

< TP(/c(r-^)) ( ^ 

< TP{Ko+Z) 



2 



Therefore, to sum up, we have 



(2.6.14) 



\9?^"'fL-m ^ PC^o) + VtP(/Co + Z) + VtF{F,H) 



Step 3: Pm case: estimate of via elliptic estimate 

For < n < A?' — 1, the n*^ order Navier-Stokes equation is as follows. 

dt{D^u^) - A^(£)t"u'") + V^(9fp™) = 



(2.6.15) 



• (D>™) = 
S"^ (£»>", aj"p'")A/' 
Dj^u™ = 



m 
in 
on 
on 



s 
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where and is given in terms of u™ and ry™ ^ . 

A straightforward apphcation of elUptic estimate reveals the fact. 

(2.6.16) < II ww'^iIl^h— n-i + iiF"iii.H— -1 + ml2H^^-^n-^;. 

A more reasonable form is as follows. 
(2.6.17) 

ll^r"'"llL2H2N-2n + l < ||^'"||2,2^2]V-2r.-l + || -ff " 1 1 ^2 Jj2N-2n- 1/2 + 1 1 ^"^ ^ 1 1 ^2 ^2]V -2n- 1 

+ wdtiD^u^ - arw")iii2^2«-2„-i + ii^r^^" - a^^^"|li2^2«-2„+l 

Then we give a detailed estimate for each term on RHS. These estimates can be easily obtained as what 
we did before, so we omit the details here. It is noticeable that the appearance of T is due to lemma 
A.15. 



|i^"||i2^2«-2„-l < TP{IC{^"^-'))ij2h 



N-2 



3=0 



2 



N-2 
j=0 



j;,oofl-2N-23-2 



n{F,H) 



< TP{Ko + Z)+H{F,H) 



N-2 

l^2^2«-2„-v.(E) ^ TP{K(r,^-^))(Y,\diu^ 



2 



N-2 
j=0 



J,oofl-2N-2j-2 



+ n{F,H) 



< TP{ICo + Z)+n{F,H) 



N-l 



< TP{)Co+Z) 

, N-2 

iii)r«'"-5r«'"iii2^2™i < TP(/c(,y™-^))^b^™ 

< TP{Kq + Z) 

Therefore, to sum up, we have 

N-l 

(2.6.18) l|5r«'"llL=H2«-2n+i < TP(/Co + Z)+ H{F, H) 

n=0 



/^oo^2iV-2j-l 



Step 4: Fm case: synthesis of estimate for 

Combining (2.6.14) and (2.6.18) with lemma A. 14 implies 

(2.6.19) Q(u™) < P(/Co) + VtP{1 +lCo+Z) + VtJ^{F, H) + n{F, H) 
Then, by forcing estimate (2.39), we have 

TiF^H'') < P{K{rr-^))+P{Q{u'^-^))<P{Z) 
niF°,H^) < T(^P(/C(r7™-^)) + P(Q(u"-^))^ <TP(Z) 

Hence, to sum up, we achieve the estimate 

(2.6.20) Q(«'") < P(/Co) + VtP{ICo + Z) 
This is actually 

(2.6.21) Q{vr) < CP(/Co) + VtCP(/Co + Z) 

for some universal constant C > 0. So we can take Z > 2CP(/Co). If T is sufficiently small depending 
on Z, we can bound Q(m™) < 2CP(/Co) < Z. 
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Step 5: P„i case: estimate of 77™ via transport estimate 
Employing transport estimate in lemma 2.37 

(2.6.22) /C(77") < P(/Co) + P(QK")) 

for P(0) = 0. Hence, we can take Z > P(/Co) + P{2CP{JCo j), then we have 

(2.6.23) /C(ry'") < Z 



Step 6: Pm case: estimate of J"^{t) 

Define ^"^ = r?™ — 770, then transport estimate in lemma 2.38 implies ||^™||^oojif5/2 ^ TZ. Naturally, 
suptg[o,T] {^"'{t) - ^^"(0)1 < llC^lliooj^s/a < TZ. Thus if we take T < 5/{2Z), then J™(t) > 5/2. A 
similar argument can justify the closeness assumption of 77™ and r}o in studying linear Navier-Stokes 
equation. 

Synthesis: 

Above estimates reveals that if we take Z — CP{K,q) where the polynomial can be given explicitly by 
summarizing all above and T small enough depending on Z, we have 

(2.6.24) S(u™)+/C(77") <Z 
and 

(2.6.25) r\t)>5/2 for te[0,r] 

Therefore, is verified. By induction, we conclude that P„ is valid for any n S N. It is noticeable that 
Z = P(/Co) where P(-) actually satisfies P(0) = 0. 

□ 

Theorem 2.44. Assume exactly the same condition as Theorem (2.43), then actually we have the 

estimate 

(2.6.26) /C(u-,p'") + JCinn < P(/Co) 
for a polynomial satisfying P(0) = 0. 

Proof This is a natural corollary of above theorem. Since we know Q(u'"-i) < Z, by remark 2.42, it 
implies IC{u"^,p"^) is bounded for any to G N. For convenience, in the following we also call this bound 
Z. □ 

2.6.4. Contraction Theorem. Define 

(2 6 27) ^(^'95^) = Ml'^h^ + 11^1112^3^+ WdivWl^^o + \\dtv\\l2Hi + Ikllioo^i + hWl^H^ 

Theorem 2.45. For j = 1,2, Suppose that , , and achieve the same initial condition for 
different j and satisfy 

dtv^ — A_4j + V ^3 = dtCbK^ dsw^ — ■ V ^3 in Cl 
V_4j ■ = in SI 

SA:i (q^ , = CW^' on S 



(2.6.28) 



v^ = on Sfe 

dtC^ = ■ Af^ in CI 



where , Af^ and are in terms of C,^ . Suppose lC{w^ ,0), K{v^ ,q^) and )C{C^) is bounded by Z. Then 
there exists < T < 1 such that for any <T <T, we have the following contraction relation 

(2.6.29) OT(w^ - v"^, -q^;T)< ]^^{w^ - 0; T) 

(2.6.30) 9Jl(C^ -(^■,T)< m{w^ - w^ 0; T) 
Proof. We divide this proof into several steps. 
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Step 1: Lower order equations 
define v 



V — V , q = q — q , w = w — w 



and C = ~ C^) which has trivial initial condition. Then 



they satisfy the equation as follows. 
(2.6.31) 



+ V^i • 5^1 ((?, I)) = -f/"^ + V^i • (D(^i_^2)z;2) in fl 

- v^ H"^ in Q. 

Sa, («, v)N^ =H^+ D(^i_^2)i;Wi on E 

V = on Eft 



where 
(2.6.32) 



V(^i_^2)(P^2i;2) „ (^1 _ ^2^Vp2 + dtC^bK^dsw + dtCbK^dzw^ 
-V(^i_^2) • -y^ 

-g2(AAi - Af^) + ©^iz;2(Ari - AA2) - D(^i_^2)t;2(AAi - J^^) + (J^^ + (^{J^^ - AA^) 
The solutions are sufficiently regular for us to differentiate in time, which is the following equation. 

(2.6.33) 



dtidtv) + V^i • .S^i {dtq, dtv) = H'^ + V^i • (D(aj^i_gj^2)u2) in Vt 

V^i • dtV = H'^ in n 

SAA9tq,dtv)Af^ =H^+B^e^A^_Q^A-)V^^^ on S 

dtV = on S(, 



where 
(2.6.34) 

= dtH^ + Va^i • (D(^i-^2)w2) + V^i • {B(^A^_A2)dtv^) + Va,A^ ■ (D^it;) + V^i(Da,^iu) - Vg^A^q 

H^= dtH^-Vo^A^-v 

= dtH^ + B^A^_A2)dtV^Af^ + D(^i_^2)v25tA/'i - {q, v)dtAf^ + Bg^A^vAf^ 



Step 2: Energy evolution for dtV 

Multiply J^dtV on both sides to get the natural energy structure: 



I [ \dtv\^J' + l f I pA^dtv\^J' = f I J'H'-dtV+ f I j'H^dtq- f I H'-dtV 
^ Jn ^ Jo Jn Jo Jn Jo Jn Jo Jt, 

+ \ f I \dtv\'dtJ^- \ f f J'B^9,A-_g^A^)v'' :l}A-dtV 
^ Jo Jn ^ Jo Jn 



LHS is simply the energy and dissipation term, so we focus on estimate of RHS. 



Jo Jn 



dtV < Z 

< VfzV^ 



f 




Jo 





HO 



\\dtv\\^o<2\\dtv\\^^^oVT 



J'H'dtq < I J'qH 
< 



l„u2 

n Jo Jn 

2 , luii 



[ jW + lklLooffo 

Jn 



s/TZ 



[ [ H^'-dtv < 
Jo Jt, 



< [ J^qH'^ + VfzV^ 
Jn 



H-i/2(s) 



iooJf-l/2(E) 



< \/f\/W 

ffldtvfdtJ' < \\dtJ'\\^^jj,T\\dtv\\l^^o 

Jo Jn 



< TZW 



WELLPOSEDNESS AND DECAYING PROPERTY OF 



VISCOUS SURFACE WAVE 



33 



So wc need several estimate on . For the following terms, it suffices to show they are bouridcid. We 
use the usual way to estimate these quadratic terms. Since we have repeatedly used this method, we will 
not give the details here. 



< z 



m\\ 



+ 11^111,2^2 + \\q\\L2m 
< Z ( Tot" + yOT» - 



^2^1 + ||w||i2^1 

< z[\f^^\fm 



i«H-l/2(S) 



^ ^( ||CIIl=ohV2 + ||5tC|LooHl/2 + IklliocHl + II^^IIl=>oH2 ^ 

< zi^\f^^\f^ 

Also we have the following estimate. 

/ J\E'' < Z\\q\\^^^o 
Jq 



L'=°H<> 



C satisfies the equation 



^ ^ IklliooHO ||Cllz,ooHl/2 + ||5tC||i=oHl/2 + ll^^lliooHl ^ 
< ZV^(^ IICILocHl/2 + l|5tCILooHl/2 + Wvh^H^ 



dtC + wldiC + wld2C = -J^^ 



. C(o) = 

Employing transport estimate and the boundedness of higher order norms, we have 

IICILocHa/2 < Vtz \\w\\^,hw, < VtzV^ 

dtC satisfies the equation 

dtidtC) + wldiidtC) + w^2d2{dtO = -M^ ■ dtw - dtN^ ■ w - dtwld^c - a^^aaC 

dtm = 
Similar argument as above shows that 

mWL^m/^ ^ VTz(\\dtw\\^,H,/. + 1^1^2^1/2 ) < VfzV^ 



Combining above transport estimate and lemma A.f4. we have the final version 

/ Jig^2 < ^zV^V^ + ZV^ \\v\\i^^jj, 
Jn 



< VtzVW'V^ + zVm" 



\\dtv\\ 



< s/fZy/^y/^+y/fZy/^(^\v\\^^H- + \\dtv\\LooHO^ 

< \/tz[VwVw' + 

Then we consider simplify the last term in RHS of energy structure. 
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where e should be determined in order for the second term in RHS to be absorbed in LHS. 
To summarize, using Cauchy inequality 

(2.6.35) \\dtv\\iooH0 + \\dtv\\l2H^ < v^^(ai"' + 0^" + 9}t) 



Step 3: Elliptic estimate for v 

based on standard elliptic regularity theory, we have the estimate 

|2 2 ^ II o Il2 , II Trll|2 , II rj-2l|2 



set r = and take L°° on both sides 

+ ||V^i • {^iA^-A-)V^)\\looHO + PiA^-A-)V^^f^looHhiJ:) 

set r = 1 and take on both sides 

\\V\\12H3 + \\q\\l2H2 < 119*^11^2^1 + \\H^\\l2Hl + \\H^\\l2H2 + 1 1 ^2 f (s) 

+ ||V^i • {n^A^-A-)V^)\\l2Hr + ||»Mi-.A=)^''-A^il2Hi(E) 

We need to estimate all the RHS terms. We can employ the usual way to estimate quadratic terms in 
L^H^ norm, however, for L°°H^ norm, we use lemma A. 14. Then we have the estimate 



RHS<\\dtv\\[2H^ + \\dtv\\[^HO+Tz(^^ +dK 



To summarize, we have 
(2.6.36) 



qWl^H- + \\v\\l^.H- + \\q\\l2H2 < 119*^1^2^. + WOtvWl^^o + TZ 01" + 



In total of (2.6.35) and (2.6.36), we get the succinct form of estimate 
(2.6.37) 0^" < y/TZ{m^ + Ot" + 9Jt) 



Step 4: Transport estimate for C, 
C satisfies the equation 

dtC + w\diC + wld2C = ■ w 



C(o) = 

A straightforward application of transport estimate can show 



■2 - will 



< VTZ\\w\\l2HW2^j:)<Vfzm^ 

Similar to the proof of transport estimate in theorem 2.37, we can use the transport equation to estimate 
the higher order derivatives. 

l|9tCllL=oH3/2(£;) < ||C^||i=ojy5/2(J^) ||^i^|lL~iJ3/2(s) + ||C|IlooH5/2(E) l|w'l|li~H3/2(s) < ZVT 

In the same fashion, we can easily show 

To sum up, we have the estimate 

(2.6.38) ^Di<'yr 
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Synthesis: 

In (2.6.37), for T sufficiently small, we can easily absorbed all D^" term from RHS to LHS and replace 
all m with m'^ to achieve 

(2.6.39) < VTzm."^ 

Certainly, the smallness of T can guarantee the contraction. 

□ 

2.6.5. Local Wellposedness Theorem. Now we can combine theorem 2.44 and 2.45 to produce a unique 
strong solution to the equation system (1.3.8). 

Theorem 2.46. Assume that (w^jvy") satisfies /Co < oo and that the initial data are constructed to 

satisfy the N^'^ compatible condition. Then there exist < To < 1 .such that (fO < T < To, then there 
exists a solution triple {u,p,r]) to the Navier-Stokes equation on the time interval [0,T] that achieves the 
initial data and satisfies 



(2.6.40) JC{u,p) + JC{t]) < CP(JCo) 

for a universal constant C > and a polynomial P(-) satisfying -P(O) = 0. The solution is unique among 
functions that achieve the initial data and satisfy IC{u,p) + /C(ry) < oo. Moreover, ri is such that the 
mapping ^ is a (72JV-2 ^iffgomorphism for each t € [0, T] . 

Proof. Since we have proved the boundedness and contraction of the iteration sequence as [1] did, then a 
similar argument as in theorem 6.3 in [1] can justify our result. Hence, we will omit the details here. The 
basic idea is that using boundedness to derive weak convergence and using the interpolation between 
lower order and higher order and the contraction to derive strong convergence such that we could pass 
to limit in the iteration sequence. □ 



3. Global Wellposedness for Horizontally Infinite Domain 

3.1. Preliminaries. In this section, we will present two types of formulation for the system and describe 
the corresponding energy structure. Since they are almost identical to that in section 2 of [2], we will 
not give detailed proofs. 



3.1.1. Geometric Structure Form. Suppose that 77, u are known and A, J\f, etc. are given in terms of 77 
as usual. Then we consider the linear equation for {v,q,^) 



(3.1.1) 



' dtV-dtfibKd3V + u-S7^v + S7^-SA{q,v)=F^ in Cl 

• u = F^ in n 

SA{q,v)Af = ^Af + F^ on S 

dt^-vJ\f = F^ on S 

V = on Sb 



We have the following energy structure. 

Lemma 3.1. Suppose {u,p,ri) satisfies system (1.3.8) and {v,q,^) is the solution of system (3.1.1), then 
(3.1.2) 

dt(\ [ J\vf + l I [ J\^Avf= I J{vF^ + qF^)+ f -v ■ F^ + ^F^ 



The geometric structure form will be utilized to estimate the temporal derivatives. Hence, we may apply 
differential operator d" = 8^° to system (1.3.8) with the resulting equation being (3.1.1) for v = d°'u, 
q = d"p and ^ = d^rj, where 
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(3.1.3) 


= 






(3.1.4) 


Fl'' = 
I 




-^dtf]d^{bK)dsUi 






O<0<a 0<p<a 




(3.1.5) 


Pi = 




8a \ 






0<l3<a ^ 




d 1 p,^ 

yo.i.o) 


pl,3 _ 


\^ (1 ar)f^ A ■iri"'~f^ri,( A- r) n ■ A- A ■ f) ii) 








0</3<a 




(3.1.7) 


Fl'' = 


y Ca,pAjkdk{d^Aiid''-^diUj + d^Ajid''-^i 








0<fS<a 




(3.1.8) 


fI'^ = 


d<^dtf,bKdzu, 




(3.1.9) 




Ajkdk{d°'AiidiUj +d°'AjidiUi) 




(3.1.10) 




p2 ^ ^2,1 _^ p2,2 




(3.1.11) 




F^'^ = - y CaBd^Aijd"-^djUi 

/ J ^ih' '■J J '■ 








0</3<Q 




(3.1.12) 




p2a ^ -d'^AijdjUi 




(3.1.13) 


f^ 


_ p3,l _|_ ^3,2 




(3.1.14) 












O<0<a 





(3.1.15) F''^ = Yl Ca,0(d^iMjAim)d"-^dmUj + d^iMjAjm)d''-^dmU^ 

O<0<a ^ ^ 



(3.1.16) F^= Y Cc,,fsd^Dr]-d'^-h 

0</3<a 

In all above, C^,/? represents constants depending on a and (i. 



3.1.2. Perturbed Linear Form. We may also take the system (1.3.8) as the perturbation of the linear 
system with constant coefficients. 



(3.1.17) 



dtu - Aw + Vp = in Q 

V • w = in 

{pi - - 77/)e3 = G^ on S 

dtil — U3 = on E 

u = on Sb 



This form satisfies the following energy structure. 

Lemma 3.2. Suppose {u,p,T]) satisfies system (3.1.17), then 

(3.1.18) 

dt(\ I \u\' + \ [ I J\Ou\^= f u-G'+pG'+ j -u-G^ + nG^ 



Compare to the original equation (1.3.8), we can write the detailed form for the nonlinear terms. 
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(3.1.19) = G^'i + + G^'^ + G^'* + G^'S 

(3.1.20) G]'' = {6ij-Aij)djp 

(3.1.21) G]'^ = ujAjkdkUi 

(3.1.22) G]'^ = [K\l + A^ + B^)-l]d33Ui-2AKdi3Ui-2BKd23Ui 

(3.1.23) G]'^ = [-K^{l + A^ + B^)d3J + AK^{diJ + d3A) 

+BK^{d2J + dsB) - K{diA + d2B)]d3Ui 

(3.1.24) G]'^ = dtfj{l+X3/b)Kd3U3 



(3.1.25) = (1 - J){diui + d2U2) + Ad2Ui + Bd3U2 



(3.1.26) 

/ p - ri - 2{diui - AKdsUi) ^ -d2Ui - diU2 + BKdsUi + AKd3U2 

G^ = diT]l -d2Ui - diU2 + BKdaUi + AKd3U2 ) +527/1 p - rj - 2{d2U2 - BKd3U2) 

^ — 9iU3 — Kd3Ui + AKd3U3 ' ^ —d2U3 — Kd3U2 + BKdsus 

, (K-l)d3Ui+AKd3U3 

+ {K~l)d3U2 + BKd3U3 



2{K - l)d3U3 



where based on the equations we have that 



(3.1.27) 

p-r] = d\r]{-diU3 - Kd3U\ + AKd3U3) + d2'n{-d2U3 - Kd3U2 + BKd3U3) + 2Kd3U3 



(3.1.28) G^ = -Dri-u 



3.1.3. Preliminary Estimates. Here we record several preliminary lemmas in estimating. 



II 2 

Lemma 3.3. There exists a universal < 6 < 1 such that if \\t]\\h^/2(j2) ^ then 



(3.1.29) ||j_i||2^ + ||^||2^ + ||B||2^< 1 



Proof. The same as lemma 2.4 in [2]. □ 

Lemma 3.4. For i = 1,2, define Ui -.T, by 

(3.1.30) Ui{x') = J J{x',X3)ui{x',X3)dx3 

Then dtri = —d\U\ — 8202 on E. 

Proof. The same as lemma 2.5 in [2]. □ 

3.2. Definition of Energy and Dissipation. In this section we will give the definition of energy and 

dissipation with or without minimum count. Compared to those in [2] (2.46-2.55), the main difference 
here lies in the full dissipation part, where we consider at most two vertical derivatives for velocity and 
one vertical derivative for pressure. It is this key improvement that largely simplifies the whole proof. 
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3.2.1. Energy and Dissipation. We define the energy and dissipation with 2-minimum count as follows. 

2 



(3.2.1) 


^n,2 = 


II n^"" 

||-'^2 




2 




(3.2.2) 


'Dn,2 = 










(3.2.3) 


^n,2 = 


\\D^u 


n 


2 

U2n-2j 





HO 



4- II D2"nlP 



n-1 



Jj2n-2j-l 



i=i 



Jj2n-2j(s) 



(3.2.4) Vn,2 = 



r +iivD2"-i»ir 

\fj2 ^ II V-'->'2 ■f'llii-o 



n+1 



2Ti-2j+2 



i=2 



//1/2(S) 



We define the energy and dissipation without minimum count as follows, where Ixu and Ixr] denote the 
Riesz potential of u and r] as defined in (A.4.1) and (A. 4. 2). 



(3.2.5) 4 = ||7Ati||^o + ||/a»7|Iho(s) + + 

(3.2.6) Vn = ||©7aw||ho + ||-Dg"Bw " 



(3.2.7) £n 



whuwio + whvfHo^^) + E ^2„-2, + E rip 

j=0 j=0 



Jlf2n-23-l 



(3.2.8) V„ 



I f\2n—l ||2 I ||v7 fS2Ti— 1_| 

I-Dq m||^2 + II V-Dq p| 



2 



+ llD^^-^nll^ +llaD2"-2„|P 
^11-^0 '/|Ihi/2(e) ^ lr*-^o '/||jfi/2(E) 



■Ei^*^'^o 

i=2 



2ri-2i+l 



Finally, we define the energy with 1-minimum count as follows. 

2 



(3.2.9) 



li^^ll 



^^. + l|i?p||^2„-2 + 5:||a,^p 



JJ2n-2j-l 



n 



+ ||Z)ry||^2„-i(j:)+^|5/r? 



H2n-2j(E) 



Since we only consider the decaying for energy and utilize the interpolation relation between £n,2 and 
£n,i to estimate, it is not necessary to define the dissipation with 1-minimum count and the horizontal 

quantities. 

3.2.2. Other Necessary Quantities. We define some useful quantities in the estimate as follows. 

2 



(3.2.10) 



^ = l|Vu||i. + ||v2n||^^+^pn,||^2 



(3.2.11) J^2N = |h||H4N+i/2(E) 

The total energy in our estimate is defined as follows. 

(3.2.12) gN{t)= sup f2iv(r)+ V sup (1 + r)'"+^£:;v+2,m(r-) + sup ^ 

0<r<t — ; 0<r<t 0<r<t U 



^2N{r) 

(1 + r) 



Note that our definition of total energy is different from that of (2.58) in [2], since we do not include the 
integral of dissipation and concentrate on the instantaneous quantities. 



WELLPOSEDNESS AND DECAYING PROPERTY OF 



VISCOUS SURFACE WAVE 



39 



3.2.3. Some Initial Estimates. Based on remark 2.6-2.8 in [2], we naturally have the following lemmas. 
Lemma 3.5. £n,2 satisfies that 



(3.2.13) 



2^1 2 "lli/0 



Lemma 3.6. For N >A, we have £n+2,2 % ^2jv and ©jv+2,2 ^ ^2iv- 
Lemma 3.7. We have the estimate ||-fA9t^||ffO(x;) ^ ll'^lljfo ^ ^2iv- 
3.3. Interpolation Estimates. 



HO(E) 



3.3.1. Lower Interpolation Estimates. In the following, we will provide several interpolation estimates in 
the form of 



(3.3.1) 
(3.3.2) 



ll^ll' < {£N+2.2)\£2Nf- 



N+2.,2) 



\£2nY- 



where 6* G [0, 1], X is some quantity, and ||-|| is some norm. 

For brevity, we will write the interpolation power in the tables below. For example, the following part 
of the table 



(3.3.3) 

should be understand as 

(3.3.4) 
(3.3.5) 





£n+2,2 


T^N+2,2 


X 


9 


K 



< {VN+2,2n£2N) 



1-9 



1-K 



When we write £n+2.2 ^ I^N+2.2 in a tabic, it means that 9 is the same when interpolating between 
£n+2,2 and £2N and between T>n+2,2 and £2n- When we write multiple entries, it means the same 
interpolation estimates hold for each item listed. Sometimes, an interpolation index r appears, we can 

always take arbitrary r G (0, 1). 

Although most of the estimates are straightforward, some of them have more complicated structures and 
needs further computation. The most common cases are as follows. 

(1) For < 61 < K < 1 



(3.3.6) 



■'N+2,2'-'2N 



OoM ^ '^N+2,2'-'2N 



1-K. 



_ cB ci-e 
< c9 ci-e 

^ ^Ar+2,2^2Af 



c9 pK — 9 cl — n 
' ^A''+2,2^Ar+2,2^27V 

cO CK-Ocl-K ^ c6 cl-0 

' ^N+2,2^2N ^2N ~ ^N+2,2^2N 



We may apply the same fashion for £n+2,2 replaced by X'jv+2,2- 
(2) For 6I1 + 6»2 > 1 



(3.3.7) 



\X\ 



~ ^JV+2,2^2Ar ^Af+2,2^2Ar ~ '-'N +2,20 j^2,2 C^N 



?l — Ol 

'2N '^N 
% £n+2,2£2N ^ ^iV+2,2 

where we utilize the bound £2n < 1- 



With this notation explained, we can state some basic interpolation estimates now. 

Lemma 3.8. The following tables encodes the power in the L°° or interpolation estimates for u with 
its derivatives either in fl or S. 



(3.3.8) 



L°° 


£n+2,2 ^ T^N+2,2 


u 


1/2 


Du 


2/(2 + r) 


Vu 


1/2 


DVu 


2/(2 + r) 


V^u 


1/2 


DV- a 


2/(2 + ,■) 





£n+2.2 ^ T^N+2,2 


u 


A/(A + 2) 


Du 


(A + l)/(A + 2) 


V?i 


A/(A + 2) 


DVu 


(A + l)/(A + 2) 


V^u 


A/(A + 2) 


DV-u 


(A + l)/(A + 2) 
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The following tables encodes the power in the L°° or interpolation estimates for ri and r] with their 
derivatives either in Q or E. 



(3.3.9) 





£n+2,2 


'Dn+2,2 


V,V 


(A + l)/(A + 2) 


(A + l)/(A + 3) 


D?7, V77 


(A + 2)/(A + 2 + r) 


(A + 2)/(A + 3) 




1 


(A + 3)/(A + 3 + r) 


diij.Oi ij 


1 


2/(2 + /■) 





£n+2,2 


'Dn+2,2 


V,V 


A/(A + 2) 


A/(A + 3) 


Drj, Vfj 


(A + l)/(A + 2) 


(A + l)/(A + 3) 




1 


(A + 2)/(A + 3) 


dtV, dtV 


1 


1/2 



The following tables encodes the power in the L°° or interpolation estimates for p with its derivatives 
either inCl or S. 



(3.3.10) 





£n+2,2 


T^N+2.2 


Vp 


1/2 


1/2 


Dp 


2/(2 + r) 


2/3 


DVp 


2/(2 + r) 


2/(2 + r) 


D'^p 


1 


2/(2 + r) 





£n+2,2 


^5^+2, 2 


Vp 








Dp 


1/2 


1/3 


D\/p 


1/2 


1/2 


D'^p 


1 


2/3 



The following tables encodes the power in the or interpolation estimates for G* . 





£n+2,2 


2^JV+2,2 




1 


(3A + 5)/(2A + 6) 


DC 


1 


1 




1 


1 


DG^ 


1 


1 


VG^ 


1 


1 




1 


(3A + 5)/(2A + 6) 


£)G^ 


1 


1 


G* 


1 


1 



(3.3.11) 





£n+2,2 


'Dn+2,2 


Gi 


(3A + 2)/(2A + 4)) 


(3A + 3)/(2A + 6) 


DG' 


1 


(3A + 5)/(2A + 6) 


G' 


1 


TO2 


DG^ 


1 


1 


VG^ 


1 


TO2 


G' 


(2A+l)/(A + 2) 


nil 


DG^ 


1 


m2 


D'^G-^ 


1 


1 


Q4 


1 


m2 



for 



nil 
ni2 



mm{{2X^ + 6A + 2)/(A^ + 5A + 6), (3A + 3)/(2A + 6)} 
: min{(2A^ + 7A + 4)/(A2 + 5A + 6), (3A + 5)/(2A + 6)} 



Proof. The estimate follows directly from the Sobolev embedding and lemma A.11-A.13, using the bounds 
II-^A^IIho ~ and ||/A9t?7||^o ^ £2N which is a natural corollary of lemma 3.7. For the terms of G*, 
we may utilize the estimating method for complicated structures stated at the beginning of this section 
and the natural product rule 



(3.3.12) \\XY\\ 

(3.3.13) \\XY 

(3.3.14) \\XY 
where in bound, we may take the larger value of 9 produced by the two options. 



2 



< 
< 
< 



\x 

\x 
\x 



2 
|2 

1 1,2 



□ 
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Then we can show several important lemmas which greatly improve the estimates above and will play 
the key role in the a priori estimates. 

Lemma 3.9. We have the estimate 

(3.3.15) k: < s^ij'^'-^s'Ji'r' 

Proof. 

2 

(3.3.16) /C = \\Vu\\l^ + \\V'u\\l^ + ^ WDuiWl,^^^ 

We need to estimate each term above. 
1. llV^lli^:^ ^ 
By the definition, we have 



1. llV^lli^:^ ^ 

finition, we have 

(3.3.17) \\yu\\l^<\iDu\\l^ + \{dsuf, 
By the divergence equation, we have 

(3.3.18) \\d3Us\\l^<\\Du\\l^ + \\Gm^ 
For i = 1, 2, by Poincare we can naturally estimate 

(3.3.19) liaauJi- < piu^l^ + ||53Wi||i^(s) 
The upper boundary condition implies that 

(3.3.20) ||a3Wi||io.(s) < P«3|li»(s) + < WDVusWl^ + \\G'\\1^^^^ 
Considering the Navicr-Stokcs equation, wc have 

(3.3.21) 1193^11'. < WdmWl^ + \\D'u4l„ + W^pWI- + PXoo 

Therefore, to summarize all above 
(3.3.22) 

llV^lli.. < \\Du\\l^ + \\Gm^ + \\DVu,\\l^ + ||Gil^(s) + \\d,u4l^ + \\D'u4l^ 
+ \\G'\\l^ + \\Dp\\l^ 
The estimates of each term above with lemma 3.8 satisfy our requirement 

By the definition, we have 

(3.3.23) \\^'4l^ < WDVuWl^ + plu^Wl^ + piusWl^ 

The first term naturally satisfies our requirement while the second term has been successfully estimated 
in the first step above. Hence, we only need to estimate the third one. Taking ^3 in the divergence 
equation reveals 

(3.3.24) \\dly4l^ < WDVu^Wl^ + \\VG%^ 
where all the terms can be easily estimated via lemma 3.8. 

3. ELi Il-Dwi||ff2(s): 

By trace theorem and Poincare lemma, we have 

2 

(3.3.25) ^ \\D''^^4lo + ||^'Vw,||5^„ + WDuiWl, 

The first and second term are naturally estimated, so we only need to focus on the last one. By Poincare, 
we have 

(3.3.26) W^Uifm % W^dz-^ifm + ll^'^^ll^o 
and also 

(3.3.27) \\Ddo,UifH. < \\Ddlui\\l, + 11-0^3^^11^0(2) 
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The boundary condition shows that 

(3.3.28) WDdsUiWlo^^^ < ||£'V||io(s) + \\DG%o^j:^ ^ ll^'^3«3||J,o + \\DG%,^^^ 
By Naiver-Stokes equation, we have 

(3.3.29) \\Ddiui\\lj, < WDdtUiWlo + \\D\i\\lo + \\D^p\\Io + \\DG'\\1, 
Hence, in total, we have 

(3.3.30) 

2 

^ \\Du,\\l, < \\D^Vui\\l, + + \\D^d3U3\\lo + WDG^'WIo^^^ + \\Ddtu4lo 

1=1 

+ \\Dml, + \\DG'\\l, 

The estimate of each term above with lemma 3.8 satisfies our requirement, so we finish this proof. □ 
Lemma 3.10. We have the estimate 

(3.3.31) £;v+2,2 < Sll^^'^'^V^.'t^r^'' 

Proof. The dissipation 'Dm+2,2 can control all the terms in f jv+2,2 except that dtt], D'^r] and D^'^+^^tj. 
A direct application of lemma A. 10 reveals that 

(3.3.32) ||z?2(iV+2)^||^^ < f V(4iV-7)2;,(4r.-8)/(4iV-7) 

Since N >3, above estimate satisfies our requirement. Also lemma 3.8 reveals that the estimate of D^r] 
will not be an obstacle. Hence, the dominating term is dtrj. By transport equation, we have 

(3.3.33) l|at7?||H0(E) ^ II«3|IhO(E) + l|Gi?,o(s) 

By Poincare, we can estimate 

(3.3.34) lk3|lH0(E) ^ WdsUaWlo 
while the divergence equation and Poincare imply that 

(3.3.35) Wd^uaWlo < \\Dui\\lo + \\G%„ < WDdsmWlo + \\G%, 
Similar to the third step in the proof of lemma 3.9, by Poincare, we have 

(3.3.36) \\Dd3Ui\\lo < \\Ddiui\\l, + \\DdsUi\\lo^^^ 
The boundary condition shows that 

(3.3.37) WDdsUiWlo^^^ < ||£>V||io(s) + \\DG%o^^^ < \\D^dsu4l, + \\DG%,^^^ 
By Naiver-Stokes equation, we have 

(3.3.38) \\Ddiui\\%, < WDdtUiWlo + \\D\i\\]jo + \\D^p\\ho + \\DG'\\]j, 
where D'^p dominates. Utilizing Poincare, we have 

(3.3.39) \\D^p\\ho ^ \\D^93p\\lo + ||^'p||io(s) 
in which we may again use the boundary condition to estimate 

(3.3.40) \\D'p\\lo^^) < ||i5\||Lo(s) + \\D'dsUs\\l„^^^ + \\D'G%,^^^ 
Hence, in total, we can estimate that 

(3.3.41) 

WdtvWlo < \\D^dsU3\\l, + ll^'G^II^,,^^^ + WDdtUiWlo + \\D\i\\lo + \\DG%o 

+ ll-D ^3^11^0 + \\D^v\\hO(^s) + ll^^^3M3||^o(E) + \\G^\\ho + \\G ||^o(s) 
+ \\D'^G^\\^„^^^ 

We may easily check via lemma 3.8 that each term above satisfies our requirement, so we finish this 
proof. □ 
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Lemma 3.11. We have the estimate 

Proof. Wc only need to estimate the lower order terms which merely appears in £n+2.i, which are Du, 
Dp and Drj. Based on the lemma 3.8, the estimates of Du and Drj have satisfied our requirement, so we 
focus on Dp. By a similar argument as lemma 3.10, we have that 

(3.3.43) \\Dp\\lo < \\Dp\\ \\mp\\ 

^ II^^IIho(s) + \\Ddiu\\l„ + \\G'\\l„^^^ + \\DAu\\l„ + \\DG'\\l„ 

where each term can be estimated directly through lemma 3.8. Then our result easily follows. □ 

3.3.2. Higher Interpolation Estimates. Then we need to give several estimates for the higher order terms. 

Lemma 3.12. We have the estimate 

f3 3 44) IIZ?2W+4 ||2 I ||v2W+5-||2 < p2/(4JV-7)25(4W-9)/(4JV-7) 

Proof. See lemma 3.18 in [2]. □ 
Lemma 3.13. We have the estimate 

(3.3.45) + ||V«||^.(j,) < 
Proof. Trace theorem implies that 

(3.3.46) \\u\\l.^^^ + \\S7u\\l,^^^ < ||V«||^o + HV^^H^o + ||^'V«||^o + H^'V^^H^o 

The last two terms naturally satisfy our requirement, so we concentrate on the first two. Similar to the 

argument in lemma 3.9, we have 

(3.3.47) llV^^ll^o + IIV^^^II^o < \\Du\\lo + \\G%, + \\DS/u\\lo + \\G%„^^^ 

+ \\dtu\\lo + \\D^u\\l, + \\Dp\\lo + \\G'\\l, + \\WG%, 
The estimate of each term on the right hand side easily lead to our result. □ 

Lemma 3.14. Let P = P{K, Drj) be a polynomial of K and Drj. Then there exists a 9 > such that 

(3.3.48) WiD^^+'^rjml,,,^^^ + \\{D^''+'n)PVu\\l,,,^^^ < eWDN+2,2 
Let Q = Q{K, b, V^) be a polynomial. Then there exists a 9 > 0, such that 

(3.3.49) \\{V^''+''fj)QVu\\l, <£l^Vr,+2,2 

Proof. This lemma is identical to lemma 3.19 in [2]. The above two lemmas implies that all the results 
used in proving lemma 3.19 have been recovered now, so this lemma is still valid. □ 

3.4. Nonlinear Estimates. 

3.4.1. Estimates of Perturbation Terms. Now we will give several estimates for the nonlinear terms in 
perturbed linear form. Similar to theorem 4.1 and 4.2 in [2], the estimating methods for these terms are 
quite standard, so we will merely give the basic rules in estimating: 

(1) Since all the terms are quadratic or of higher order, then we may expand the differential operators 
with Leibniz rule and estimate resulting quadratic terms by the basic interpolation estimate in 
lemma 3.8 either in or L^ norm, combining with the definition of energy and dissipation, 
trace theorem and Sobolev embedding. For the choices of form (3.3.13) or (3.3.14), we should 
always take the higher interpolation index for En +2.2 and 'Dn+2.2- 

(2) For the 2N level, there is no minimum count for the derivatives, hence we only need to refer to the 
definition of energy and dissipation; however, m N + 2 level, we have to resort to lemma 3.8 for 
terms that cannot be estimated directly. Also, lemma 3.14 can be utilized for more complicated 
terms. 

(3) Note that the most important difference between our proof and that of theorem 4.1 and 4.2 in [2] 

is that we concentrate on the horizontal derivatives instead of the full derivatives. Thereafter, in 
the estimates related to X'jv+2,2, we should try to avoid the introduction of vertical derivatives 
via Sobolev embedding, especially for pressure p. 
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(4) For the term in form of !l-^|l//2(2)> instead of using trace theorem directly, we should first write 
it into 11-^^11^0(5]) + ||-DX||^o(2) + 1 1 ■^^"'^11 HO (s) then apply trace theorem, which will only 
produce one vertical derivatives rather than three. 

Considering that the proofs for the following estimates are lengthy but not difficult based on the above 
rules, and also most of them have been achieved in that of theorem 4.1 and 4.2 of [2], we will omit the 
details here and only present the theorems. 

Lemma 3.15. There exists a 9 > such that 
(3.4.1) 



^2^2(A,+2)-4^1 



HO 



+ 



2yj2{N+2)-3q2 







HO 



+ 



^2^2(iV+2)-4g3 



and 
(3.4.2) 

D. 



2(JV+2)-l^l 



HO 



^2(iV+2)-1^2 



//I 



D 



2(jV+2)-1^3 



Hi/2(E) 



^2(iV+2)-1^4 



< S2nT^N+2,2 



Lemma 3.16. There exists a 9 > such that 
(3.4.3) 



|V74Ar-2^1||2 ||v74JV-1^2||2 1 1 y74W-2^3 1 1 ^ ^ cl+0 



and 
(3.4.4) 



^4W-lg2| 



\ho 11-^0 



I _L II n4W-1^4||^ 

l//i/2(E) 11^0 ^ llffi/2(s) 



3.4.2. Estimates of Other Nonlinearities. Here is a few lemmas for some other nonlinear terms which 
have been proved in proposition 4.3-4.5 of [2]. Although our definition of energy and dissipation is a 
little bit different, we do not even need to revise the proof here. 



Lemma 3.17. We have that 

|2 



(3.4.5) 
and 
(3.4.6) 
Also 
(3.4.7) 

Lemma 3.18. It holds that 



IxGXh^ + WIxGXh- + WhdtGXjj, < S2Nrmn{e2N,V2N} 
ll-fAG^IIjji(s) + II-^aG'^IIjji(s) ^ ^2Jvniin{f2JV,^'2Jv} 



(3.4.8) 



and 
(3.4.9) 
Also 
(3.4.10) 

Lemma 3.19. We have that 
(3.4.11) 

and 

(3.4.12) 



\\h[{AK)d3U^ + {BK)d:iU2]\\l, + J2 \\h[udiK]\\lo < V\ 



2 

2JV 



\h[{l-K)G-']f^,<£2NVl^ 



\dr^'Af^,<v2M 



|5f+^^||^„<P;V+2,2 



WELLPOSEDNESS AND DECAYING PROPERTY OF 



VISCOUS SURFACE WAVE 



45 



3.5. Energy Estimates. In tliis section, we will present the energy estimate for horizontal derivatives. 
Since this part has been carefully completed in section 5 of [2] and is not the main improvement of our 
new proof, it is not necessary for us to show all the details. Hence, we will not give all the details of the 
proofs and only comment when we need different technique to achieve the estimates due to the changing 

of definition of energy and dissipation. 

3.5.1. Estimates of Highest and Lowest Temporal Derivatives. 

Lemma 3.20. Let 9" = ^^JV pt defined as in geometric structure form. Then 

(3.5.1) \\Fm, + \\dt{JF')\\l, + \\F' 

Proof. This is identical to theorem 5.1 in [2], so we omit the proof here. 
Lemma 3.21. Let d" = d^~^'^ and let he defined as in geometric structure form. Then 

(3.5.2) \\F^f^, + \\dt{JF^)f^, + \\F%, 
Also, if N >3, then there exists a 6 > such that 



□ 



+ 11^^ 



(3.5.3) 



\H0 ~ ^2N^N+2,2 

Proof. This is a revised version of theorem 5.2 in [2]. The first part is exactly the same as the original 
theorem, so it is naturally valid. However, we will utilize different techniques in the second part. By the 
definition, F^ = F^'^ + F^'^. By the same argument as (5.5) in [2], we have 

;i2,l ||2 



I H° ~ ^2N^N+2,2 



(3.5.4) 

For F^'^, a calculation reveals that 

(3.5.5) ^2,2 ^ d^+^(^d^f^lK)d3Ui + d^+^{d2fjbK)d3U2 - d^+^KdsUs 

We can estimate 



(3.5.6) 
Lemma 3.9 implies that 
(3.5.7) 



< 



HO 



\d3Ui\ 



< llVull 



< c'V(2+r)o2/(2+r) 
^ <--2W <--7V+2,2 



Also, we have that for < \a\ < N + 2, there exists a. 9 > such that 

(3.5.8) Wdrd^mlo < \\dMlm(^) < WdMnoi^) WdMmii:) < 4/^^+2,: 

Hence we have 



(3.5.9) 



di'+'id^f,bK) ^^^^ < ai^ii^o 



N+2-a 



ibK) 



^ ci-e cB 

~ ^2JV ^JV+2,2 



Then when r is sufficiently small, we can naturally estimate 

2 



d^+\difjbK)d3Ui 



d^+\d2f]bK)d3U2 



(3.5.10) 

Similarly, we have 
(3.5.11) 

The same argument as (5.11) in [2] shows that 

(3.5.12) Wdt^'+'KOau^Wlo 
Therefore, our result easily follows. 



HO 



< £2N^N+2,2 



< £2N^N+2,2 



HO 



S. £oiv£i 



2N^N+2,2 



□ 



The following three lemmas is just the restatement of proposition 5.3-5.5 in [2], then we will only present 

the lemmas without the proofs. 



Lemma 3.22. There exists a 9 > such that 
(3.5.13) 



\dru 



2 

HO 



+ \\drvm 



H°{-E) 



+ 



f 

Jo 



\m^''u\\^,<e2N{o) + {e2N{t)) 



3/2 



+ 



Jo 



^2JV^2JV 
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Lemma 3.23. Let be given with df = 9^"*"^. Then it holds that 
(3.5.14) 

Lemma 3.24. It holds that 

(3.5.15) \Ht)\\lo + 



l_f/0(S) 

3.5.2. Estimates of Mixed Derivatives. 

Lemma 3.25. Let a e fee 5mc/i that \a\ = 4N. Then 



^0 



(3.5.16) 



i 



E 



^ V ^2n'T>2N + \JT)2N^^2N 



Proof. Essentially, this lemma is completely identical to lemma 6.1 in [2]. We only need to notice that 
although our definition of energy and dissipation is different from that in [2] , the terms related to ry is 
exactly the same. Also we should refer to the general rules for estimating in the section of nonlinear 
estimates to deal with ||-'i^||/f2(5]) form. Hence, this lemma is still valid now. □ 

Lemma 3.26. Suppose that a € N^"*"^ is such that < 2N — 1 and 1 < |q;| < 4iV. Then there exists a 
6 > such that 

(3.5.17) \\D-^ u\\^o + \\DD u||^o + ||-Di rjWjjo^j^s^ + \\DD v\\hocs) 

+ / \\Dt^-'m\\l, + \\DD^^-'Bu\\l, < f2w(0) + / £^nV2n + VT^2nICT2n 
Jo Jo 

Proof. This lemma is identical to proposition 6.2 in [2], and the proof is almost the same. We only need 
to notice that now we should utilize lemma 3.16 to estimate the nonlinear terms to replace lemma 4.2 
in [2]. Hence, the result still holds. □ 

Lemma 3.27. We have the estimate 

(3.5.18) < 
Also, there exists a 9 > such that 

Coriq-v Iln2^4||2 ^ c'O ^l+l/(A+3) 

l-J.O.J-yj \\J^ ^ \\hO(i:)^^2N^N+2,2 

and 

(3.5.20) \\D'G^l<£|^v];^i^,'+'^ 

Proof. This is basically lemma 6.3 in [2], the only difference is that we need to estimate ||Vu||jjo directly 
now. By a similar argument as in lemma 3.9, we have 

(3.5.21) 

llVwIlLo < \\Du\\lo + \\G%„ + WDVmWlo + WG'WIo^j:^ + WdmWlo + ||l?'«i||^o 
+ \\G'\\lo + \\Dp\\lo 

Hence, an application of lemma 3.8 implies that ||Vu||^o ^ X'j^^2 2''^^^^^. Using exactly the same 
argument for the other part of the proof, we can see the validity of this lemma. □ 

Lemma 3.28. Suppose that a G is such that ao < N + 1 and 2 < |a| < 2{N + 2). Then there 

exists a 6 > such that 

(3.5.22) d,l^\\Dr+\\\l, + \\DD'^+\\\l, + \\Dr+^\l„^^^ + \\DD'^+ml„^^^ 

+ WDl^'+'BuWl,, + \\DD'^+'Bu\\l„ < £'2n'Dn+2.2 

Proof. This is almost the same as proposition 6.4 in [2], so we omit the detailed proof here. We only need 
to note that for higher derivatives, we should use lemma 3.15 to replace lemma 4.1 in [2] to estimate 
nonlinear terms. Then for lower derivatives, we should utilize lemma 3.8, 3.10 and 3.27 to replace 
lemma 3.1, 3.16 and 6.3 in [2], which share exactly the same results. Therefore, the conclusion naturally 
holds. □ 
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3.5.3. Estimates of Riesz Potentials. In the following, wc will record several lemmas related to Ricsz 
potentials, which have been proved in lemma 6.5-6.7 in [2]. Note that all the preliminary lemmas have 
been recovered in lemma 3.17 and 3.18, so we do not even need to revise the proofs. 



WIxPWho ^ ^2JV 



2 < ^A/(1+A)^l/(1+A) 

hphG^ 



HO ~ •^2N ^2JV 



2N'^2N 



Lemma 3.29. It holds that 
(3.5.23) 
(3.5.24) 

Lemma 3.30. It holds that 
(3.5.25) 

Lemma 3.31. It holds that 

(3.5.26) ||7;,«||^o + ||/a?7|Iho(e) + £ II^^awII^o < £2n{0) + £ ^^2^2 

3.5.4. Synthesis of Energy Estimates. Assembling lemma 3.22, 3.23, 3.24, 3.26, 3.28 and 3.31, we have 

the following energy estimate. 

Proposition 3.32. Let F"^ he given with d^^"^. Then there exists a > such that 

(3.5.27) dt (^£n+2,2 - 2 ^ Kd^+^pF'^^ + Vn+2,2 < SUT^n+2,2 
Proposition 3.33. There exists a 6 > such that 



-'2N 



(3.5.28) 



£2N{t) 



'D2N{r)di < £2n{0) + {£2N{t)) 



3/2 



[\£2N{r)y 
Jo 



2?2jvdr 



+ 



/ V^2N{r 

JO 



)K{r)T2N{r)dr 



3.6. Comparison Estimates. In this section, we will prove several comparison results between hori- 
zontal derivatives and full derivatives, both in energy and dissipation. 

3.6.1. Dissipation Comparison Estimates. 

Proposition 3.34. There exists a 6 > such that 



(3.6.1) 

Proof. We define 

(3.6.2) 3;;v+2,2 



V 



N+2,2 



N+2,2 



2N^N+2,2 



D 



2{N+2)-1q1 



HO 



D. 



2{N+2)-1^2 



+ 



^2(iV+2)-1^3 



^2(Ar-|-2)-1^4 



//i/2(E) 
2 

Hi/2(E) 



and 

(3.6.3) Z = ©JV-l-2,2 + yN+2,2 

Then the proof is divided into several steps. 
Step 1: Application of Korn's inequality. 

Since any horizontal derivatives of velocity vanishes on the lower boundary, we can directly apply Korn's 
inequality (see lemma A. 5) to achieve that 



(3.6.4) 



D. 



2{N+2) 



< 



D. 



2(N+2) 



HO 



V 



N+2,2 



< Z 



Step 2: Initial estimates of velocity and pressure. 

Let < j < (iV 2) - 1 and a G be such that 2<2j-F|a|<2(A/'-|-2)-l. Hence, we naturally have 



(3.6.5) 



HO 



< 



<z 



We apply the operator d"dfds to the divergence equation in (3.1.17) and rearrange the terms 

(3.6.6) d"did3{dsU3) = d'^did^i-diui - d2U2) + asG' 



48 



LEI WU 



d"dldlu3 



< 



HO 



2(JV+2) 



D. 



d"d{G^ 



Thus we have 
(3.6.7) 

Hence, it is easily imphed that 
(3.6.8) 

Apply the operator d°'dl to the third equation in (3.1.17), then we have 



< Z 



<z 



(3.6.9) 



< 



d^'dfAu 



HO 



+ 



HO 



<2 



So we only need to estimate the terms related to d^Ui and dip for i = 1, 2. 



Step 3: Estimates of velocity and pressure in the upper domain. 

In order to achieve above estimate, wo need to employ a localization argument. Define a cutoff function 
Xi = Xiixs) given by xi S C~(R) with xi^a) = 1 for X3 € Qi = [-26/3,0] and xiixs) = for 
0:3 ^ (-36/4,1/2). 

Define the curl of the velocity w = V x u, which means Wi = d2U3 — d^Ui and W2 = dsUi — diu^. 
Therefore, Xi^i ioi i = 1,2 satisfy the equation 

(3.6.10) 

Ad'^dUxim) = xid'^di+'w, + 2id3Xi){d'' didswi) + {dixi){d''diwi) - xiV x {d'^diG') 
in as well as the boundary condition 

r d'^dlixiwi) = 2d2d°'diu3 + d°'dlG^ -62 on E 

(3.6.11) I d'^di{xiW2) = -2did'^diu3-d'^diG^-ei on S 

( d^diixiwi) = d^di{xiW2) =Q onSb 
Based on the standard elliptic estimate for Poisson equation, we have 



(3.6.12) 



d'^diixiwi) 



< 



Ad'^diixiw, 



H- 



d^dlixiWi) 



ffi/2(E) 



Hence, we have to estimate each term on the right hand side. 

Let (f> e H^{n). For a ^ we may write a = /? + (a — /3) with = 1. Then an integration by parts 

reveals that 



(3.6.13) 



/ 

Jq 



cl>xid''di+'wi 



d^(j)Xid°'~'^d- 



Wi 



< Mh^ 



XiD: 



2{N+2) 



Wi 



HO 



Since 2{j + 1) + |a - /3| e [3, 2{N + 2)], we naturally have 



(3.6.14) 



XiD 



2(JV+2) 



< 



HO 



D 



2{N+2) 



Then if we take the supreme over all 4> such that ||(/>|lj:^i < 1, then we may get 

2 



(3.6.15) 



Xid'-dl+'w, 



< Z 



H- 



A similar argument without integration by parts shows that the above result is also true for a = since 
in this case j < {N + 2) — 1 implies 4 < 2{j + 1) < 2{N + 2). In a similar fashion, we may apply 
integration by parts for ^3 for the other terms in norm and achieve that 

(3.6.16) 



2{d3Xi){d'' dfdswi) 
(3.6.17) 



<{\\dsXi\\loo + \\dlxi\\lj 



n2(JV+2) 

^2 



Wi 



< 



HO 



<z 



H^ 



(3.6.18) 



{dlxiKd'^diw.^ 
XiV X {d"diG^) 



2 

2 



~ 11^3X1 llL 



< 



HO 



<Z 



H- 



<(IIXi|li~ + l|53Xi|li=») 



^2(^.+2)-l^l 



HO 



< z 



Above four estimates together yield 
(3.6.19) 



Ad^diixiwi) 



<z 



H- 
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Then for the boundary terms, a direct apphcation of trace theorem shows that 



(3.6.20) 



/fl/2(S) 



d"dId2U3 



< 



i/l/2(S) 



D. 



2(iV+2)^ 



< z 



(3.6.21) 

Then above imphes that 
(3.6.22) 





2 

< 







D 



2(JV+2)-1^3 



i/l/2(S) 



<2 



Hi/2(E) 



<2 



Hence, the eUiptic estimate gives us the final form 

2 



(3.6.23) 



< 



<z 



Since d^d'^dfui — d^d'^dKw^ + di_uz) and d^d'^dfu2 — d3d'^df{d2U3 — wi), we can deduce from above 
that for i = 1, 2, 

2 

HO{ni 



(3.6.24) 



<2 



which, by considering the Navier Stokes equation, further imphes 

2 



(3.6.25) 



did"dip 



Combining aU above, we have achieved that 



(3.6.26) 



D 



2(iV+2)-l 



< Z 

2(JV+2)-l 



D: 



Vp 



< Z 



Step 4: Estimates of velocity and pressure in the lower domain. 

Now we extend above estimate to the lower domain ^2 = [—b,—b/3]. Define a cutoff function X2 S (M) 
such that x^ixa) = 1 for xs G fl2 and X2{x3) = for ^ (—26, —b/6). Notice that this cutoff function 
satisfies that supp(Vx2) ■= f^i^ which will play a key role in the following. 
Then the localized variables satisfy the equation 

-A{x2u) + V(X2P) = -dt{x2u) + xiG^ in ^^ 

V • (X2U) = X2G2 _^_jj2 

X2U = on S 

X2U = on Sb 



(3.6.27) 



where 

(3.6.28) = 53X2 (pea - 2^3^) - (532x2)^ = ^3X2^3 

Let < i < (iV + 2) - 1 and a e such that 2 < 2 j + \a\ < 2{N + 2) - 1. Then we may apply 
the differential operator on (3.6.27) and consider the elliptic estimate (B.2.3) for Navier Stokes 

equation, which implies that 



(3.6.29) 



d"didi{x2u) 



+ 



d"didi{x2P) 



< 



HO 



d"di+'di{X2U) 



+ 



d"didi{x2G' + H'] 



d"didi{x2G^ + H^) 



2 



HO 



Similar to the argument in the upper domain, via integration by parts, it is easy to deduce the H ^ 
estimates that 



(3.6.30) 



2 



d'^didi{x2G': 



d^d{+'di{X2U) 

For the remaining part, we may directly estimate 
(3.6.31) 

< 



H-1 



+ 



d<^didi{x2G^ 



HO 





2 

+ 




2 

< 

HO 




2 

+ 

HO 







2 

+ 




i?2(ni) 



D 



2(Ar+2)-l 



HO{Q,i) 



2 

HO 
<Z 
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where we utilize the property for Vx2 stated above. Thus we have 



(3.6.32) 



did" dip 



< 



<2 



HO 



which, based on the Navier-Stokes equation (3.1.17), further implies 
(3.6.33) 



<Z 



Combining all above, we have 



(3.6.34) 



D 



2(iV+2)-l 



H°(Q.2) 



< z 



Step 5: Estimates of free surface. 

With above estimates in hand, we may derive the estimates for rj by employing the boundary condition 
of (3.1.17), i.e. 

(3.6.35) r? = P - 393^3 - Gl 

We may take horizontal derivative on both sides of above equation, i.e. for a e and 3 < |a| < 

2{N + 2) - 1, we have 

1 2 ^ iioa,. Ii2 I \\ cxa. c\ Ii2 , llqa/^SM^ 



(3.6.36) < l|5>|lW(E) + P"53t^3|lHV2(E) + ||5"G^Lv2(E) 



^^2(iV+2)-lp 



HO 



D 



2(iV+2)-l 



3|Ihi/2(e) 



< 

< Z 

For the temporal derivative of r], we need to resort to the transport equation on the upper boundary, i.e. 
(3.6.37) dtv = M3 + 

For a G and 1 < |a| < 2(A'' + 2) — 1, by trace theorem, Poincare theorem and the divergence equation 
V ■ u = G^, we have 



(3.6.38) 



l|5"5t,?||^V.(E) < 

< 
< 
< 
< 



d"u3\\l^+z 

d"u3\\lo + \\d"Dus\\lo + Wd'^dausWlo + Z 

d"d3U3\\l„+Z 



< z 



Similarly, for j = 2, . . . , (N + 2) + 1, we may apply Sj. to transport equation and for a € N and 
< |a| < 2(Ar + 2) - 2j + 2 estimate 



(3.6.39) 



Hi/2(E) 



< 



< 



2 

+ Z 



+ 



Hi/2(E) 



< z 



Step 6: Synthesis. 

To summarize all above, since f2 = Oi U 02 we have the estimate 

(3.6.40) 'Dn+2,2 ^ T^N+2,2 + 3^iV+2,2 

By lemma 3.15, it is obvious that there exists a ^ > such that 

(3.6.41) 3^;V+2,2 < <?2JV^'jV+2,2 

Therefore, our result naturally follows. 
Proposition 3.35. There exists a 6 > such that 

(3.6.42) V2N < f>2N + S2N^2N + ICJ^2N 



□ 
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Proof. It is easy to sec that most part of the proof is identical to that of proposition 3.34, except that 
we do not need the minimum count and interpolation now. Also we will utilize lemma 3.16 instead to 
estimate the nonlinear terms in perturbed linear form. Finally, a direct estimate for Riesz potential term 

(3.6.43) \\hu\\l^ < \\Bhu\\lo 

can close the proof. □ 
3.6.2. Energy Comparison Estimates. 
Proposition 3.36. There exists a 6 > such that 

(3.6.44) £n+2,2 ^ £n+2,2 + ^2N^N+2,2 
Proof. We define 

(3.6.45) 



2,2 



,2v72(iV+2)-4^1 







HO 



+ 



2^2{N+2)-Sq2 







+ 



,2v72(iV+2)-4^3 







and 

(3.6.46) Z = £n+2,2 + Wjv+2,2 

Then we can directly apply elliptic estimate (B.2.10) in perturbed linear form. For a e and |a| = 2, 

we have 

(3.6.47) + ||a>||^.(«+„_3 < \\d^dtu\\l,,.^,,^, + ||9"Gi ||5,2(«+2)-4 

+ ||9"G^||jj2(iV+2)-3 + II^"^IIh2(N+2)-7/2(J]') + 1 1 9"G"^ 1 1 ^2(N+2) -7/2 ( J^J 

Naturally, we have 

(3.6.48) ||9"G^ ||^2(N+2)-4 + ||f^"G^||^2(„_,.2)_3 + ||9"G"^||jj2(Ar+2)-7/2(J]) ^ WjV+2,2 ^2 

and 

(3-6.49) ||9"?7||^2{iv+2)-7/2(s) < ^Ar+2,2 < 2 

So the only remaining term is 

(3.6.50) ||l9"l9tU||^2(JV + 2)-4 < ||l9tU||^2(JV + 2)-2 

which should be estimated with a finite induction in the following. 
For j = 1, . . . , (A'' + 2) — 1, we have the elliptic estimate 

(3.6.51) 



dfu 



Jj2(,N + 2)-23 
+ 



+ 



< 



Jj2(N+2)-23-l 



Jj2{N+2)-2j-2 



+ 



Jj2(N+2)-23-2 



^2(JV + 2)-2j-l 

By the definition, it is easy to see 
(3.6.52) 



JJ2(N + 2)-2,--3/2(s) 



diG' 



//2(JV + 2)-2,--3/2(2) 



diG' 



fJ2(N+2)-2j-2 



dlG^ 



fJ2(N+2)-2j-l 



and 

(3.6.53) 

When j = (iV + 2) - 1, we have 
(3.6.54) 



dfr] 



H-2(N+2)-23-3/2(J]) 



d-'G^ 

* iJ2(iV + 2)-23-3/2(J]) 

^ £n+2,2 ^ -2 



< WjV+2,2 < 2 



\d^+\\\^,<£M+2,2<2 



fj2(N+2)-2j-2 

Hence, the elliptic estimate shows 

(3.6.55) ||af+^n||^2 + ||af+V|r^.<^ 

Inductively, wc can estimate all the temporal derivatives of velocity and pressure for j = 1, . . . ,N, which 
will finally close the proof. 
Therefore, we have shown that 



(3.6.56) 



^iV+2,2 ^ ^iV+2,2 + Wjv+2,2 
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Based on lemma 3.15, there exists a ^ > such that 

(3.6.57) 'Wn+2,2 < S2N^N+2,1 

Hence, our result naturally follows. □ 
Proposition 3.37. There exists a 9 > such that 

(3.6.58) £2N < £2N + 

Proof. It is easy to see that the proof is almost the same as that of proposition 3.36 except that there is 
no minimum count now and we should utilize lemma 3.16 instead. So we omit the proof here. □ 

3.7. A Priori Estimates. First, we concentrate on the 2A'' level a priori estimate. 

Lemma 3.38. There exists C > such that 

(3.7.1) sup J"2jv(r) < exp (c [ ^/^dr 

0<r<t V ^0 

^2Ar(0) + i^ (l+^2iv(r))p2iv(r)dr+ (^^ \/^(r)J-2jv(r)dr^ 

Proof. This lemma is almost identical to lemma 9.1 in [2], so we omit the detailed proof here and only 
outline some differences here. Based on trace theorem, we have that for < < 4 A'' and /3 e 

(3-7.2) P^Alu^^) ^P^4m ^'^^N 

Hence, although our new definition of dissipation is slightly different from the original paper, we can still 
estimate and achieve the same form. □ 

Proposition 3.39. There exists a universal constant < ^ < 1 such that if Qn{T) < S, then 
(3.7.3) sup J-2iv(r) < 7"2jv(0) + t / I?2jv(r)dr 

0<r<t Jo 

and 
(3.7.4) 



(3.7.5) 

for allO<t<T 



f /C(r)J-2iv(r)dr < 5(8+2A)/(8+4A)_^^^^o^ ^ ^(8+2A)/(8+4A) /•* V2N{r)Ar 
Jo Jo 

Jo Jo 



Proof. By lemma 3.9, we have that /C < £n+2'2^ ^2n'^~^^^ ^^'^ < 1, we may take r = 2A/(4 + A) to 
achieve that 

(3.7.6) /C<4^+2'2^/(«+^"' 

Then utilizing the decaying of fjv+2,2, the other part of the proof is identical to proposition 9.2 and 
corollary 9.3 in [2], so we will omit it here. □ 

Theorem 3.40. There exists a universal constant < 5 < 1 such that if Qn{T) < 5, then 

(3.7.7) sup £2N{r)+ f V^n + sup < £2n{0) + J^2n{0) 

0<r<t Jo 0<r<t '- + f 

for allO<t<T. 

Proof. This is identical to theorem 9.4 of [2]. Notice that in the section of energy estimates and com- 
parison estimates, we have recovered all the necessary results for 2A'' level a priori estimate here. □ 

Then we consider the a priori estimate for A'' + 2 level. 

Lemma 3.41. Let be defined with d^'^'^. There exists a universal constant < 5 < 1 such that if 
Qn{T) < 6, then 

(3.7.8) |^JV+2,2(i) < £n+2,2 - 2 / J{t)d^+h{t)F^{t) < ^^JV+2,2(t) 

for allO<t<T. 

Proof. It is the same as lemma 9.6 of [2]. □ 
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Theorem 3.42. There exists a universal constant < 5 < 1 such that if Qn{T) < 5, then 

(3.7.9) sup (1 + rf+^£N+2a{r) < <?2iv(0) 

0<r<t 

for allO<t<T. 

Proof, utilizing lemma 3.10, we may employ the same argument as theorem 9.7 in [2] to show this. 
Also notice that all the necessary lemmas have been recovered in the section of energy estimates and 
comparison estimates. □ 

Next, we will give an interpolation argument to achieve the decaying results for energy with different 
minimum count. 

Theorem 3.43. There exists a universal constant < 5 < 1 such that if Gn{T) < 6, then 

(3.7.10) sup (1 + r)'+^£N+2Ar) < £2NiO) 

0<r<t 

for allQ<t<T. 

Proof. By lemma 3.11, we have that 

Combining with theorem 3.42 and theorem 3.40, we can easily see that 

/ X (l+A)/(2+A) / ^ l/(2+A) 

(3.7.12) {l + ry+^£N+2,i{r) < Ul + rf+^£N+2,2{r)] K2jv(r) 

(l+A)/(2+A) . X l/(2+A) 



< U2Jv(0)j U2n{0)] =£2n{0) 



Hence, our result easily follows. □ 

Finally, we can sum up all above to achieve the a priori estimate. 
Theorem 3.44. There exists a universal constant < 5 < 1 such that if Qn{T) < 5, then 
(3.7.13) g^it) < £:2jv(0) + J-2Jv(0) 

for allO<t<T. 

Proof. To sum up the conclusion of theorem 3.40 , 3.42 and 3.43 easily implies this. □ 

3.8. Global Wellposedness. Since we introduce several terms related to Riesz potential in the a priori 
estimates, we have to revise our local wellposedness result here to adopt this, which is perfectly completed 
in theorem 10.7 of [2]. 

Theorem 3.45. Suppose that initial data are given satisfying the 2iV*'' compatible condition and ||uo||^4jv + 
||r?o||fl-4]v+i/2(s) + II-^aWoII^o + ||Ja??o||//o(s) < 00- Let e > 0. There exists a So = 6o{e) 

(3.8.1) To = C(e)min{l,.;; 1>0 



11% II 



JJ4]V + 1/2(J]) 



where C(e) > is a constant depending on e, such that ifO<T<To and ||uo||^4n + ||??o||h*'V(e) < ^O; 
then there exists a unique solution {u,p,r]) to system (1.3.8) on the interval [0,T] that achieves the initial 
data. The solution obeys the estimates 

(3.8.2) sup £2N{t)+ sup ||/ap(*)IIho + / T^2N{t)dt 

0<t<T 0<t<T Jo 

+ [ ( ||5r+'«(*)||(,Hi). + ||5rMi)||H0 )dt < e + WhuoWlo + 

and 

(3.8.3) sup ( £2N{t) - \\I\u\\]jo - ||/a??||ho(s) ) < e 

0<t<T \ J 

(3.8.4) sup J"2iv(t) < CT2n{Q) + e 

0<t<T 
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Finally, wc come to the conclusion for global wcllposcdness. Since we have recovered the a priori 
estimate and do not improve the argument for this part, the reader may directly refer to theorem 11.1 
and 11.2 in [2], which completes the whole proof. 

Theorem 3.46. Suppose that the initial data {uq, t]o) satisfy the 2N*^ compatible condition. There exists 

a K > such that if ||'Uollff4jv + |1»7o|1^4n+i/2(-5]') < k, then there exists a unique solution {u,p,rj) on the 
interval [0, oo) that achieves the initial data. The solution obeys the estimate 

(E)) 

Appendix A. Analytic Tools 



(3.8.5) Sjv(oo) < C(||uo||h4n + II%IIh*n+i/2(s)) < Cn 



A.l. Products in Sobolev Space. We will need some estimates of the products of functions in Sobolev 
spaces. Since these results have been proved in lemma A.l and lemma A. 2 of [1], we will present the 

statement of the lemmas here without proof. 

Lemma A.l. Let U denote either S or f2. 

(1) iei < r < si < S2 be such that si > n/2. Let f G H^'{U), g G H^^{U). Then fg G H''{U) and 

(A.1.1) WfgWnr < WIWh^. \\9\\h^. 

(2) LetO<r<si<S2 be such that S2 >r + n/2. Let f e H^^{U), g e H^^{U). Then fg G H''{U) 
and 

(A.1.2) < 11/11^., Mh.. 

(3) LetO<r<si<S2 besuchthats2 > r+n/2. Let f G g G H'^^iE). Thenfg G //"-^^(I]) 
and 

(A.1.3) WfgWH-^i < UWh-^ M\h^2 

Lemma A.2. Suppose that f G C^iT.) and g G H^/^{T.). Then fg G and 
(A.1.4) ll,/.9llffi/2 < ll/llci 

A. 2. Poincare-Type Inequality. We need several Poincare-type inequality in ft. Since all these lem- 
mas have be proved in lemma A.10-A.13 in [2], we will only give the statement here without proof. 

Lemma A. 3. It holds that 

(A.2.1) 11/11^2(0) <ll/lli2(E) + 11^3/11^2(0) 

for all f G H^ifl). Also, if f G W^'°°{n), then 

(A.2.2) ||/||i^(0) < + 1193/11^00(0) 

Lemma A. 4. It holds that ||/||//o(£) < ll^^a/llffo for f G H^{fl) such that / = on S^. It also holds 
that ||/||ioc(s) < 11^3/1^00(0) for f G W''°°{n) such that / = on Ef,. 

Lemma A. 5. It holds that \\u\\jji < ||Pu||jjo for all u G H^{Q,; B?) such that f = on S;,. 

Lemma A.6. It holds that ||/||jji < ||V/||jjo for all f G II^{Q) such that / = on S;,. Also, 
||/|lvKi,oo(o) < l|V/||ioc(o) for all f G W^^'^ift) such that / = on Hi,. 

A. 3. Poisson Integral. For a function / defined on S = R'^, the Poisson integral Vf in R'^ x (— oo,0) 
is defined by 



(A.3.1) rf{x',xs)= [ m 



g27r|5|a:3g27ria;'-{^^ 



where /(^) is the Fourier transform of f{x') on R"^. Then within the slab R^ x (—6,0), we have the 
following estimate based on lemma A. 5 in [2]. 



Lemma A. 7. The Poisson integral satisfies that for q gN, 

Ih" ~ ll/llij-9-l/2(S) 
llij<!(S) 

where denotes the usual homogeneous Sobolev space. 



(A.3.2) ||V«P/||^o<"""' 
(A.3.3) ||V''P/||^o<"^"' 
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A. 4. Riesz Potential. For a function / defined in fl, we define tlie Riesz potential 

(A.4.1) hf{x',xs)^ f" [ /U:C3)|^r^e2--'-«dedx3 

J-b Jr^ 

Similarly, for / defined on E, we set 

(A.4.2) hf{x')= [ /(C)|^r'e^--'-€d^ 

We have the following lemmas to describe the product of Riesz potential and its interaction with the 
horizontal derivatives as lemma A. 3 and A.4 in [2]. 

Lemma A.8. Let A e (0, 1). If f e H°{Q) and g,Dg e il^(O), then 

(A.4.3) \\h{f9)\\H0 < UWho \\9\\h^ WDqWh' 

Iff€ H°(S) and g,Dg e H^S), then 

(A.4.4) I|/a(/5)IIho(e) ^ II/IIho(e) II5|Ihi(e) WDgWmU 

Lemma A. 9. Let A e (0, 1). /// G H''{fl) for k > 1 an integer, then 

(A-4.5) II^Ai^VLo ^ ||^'"V||io \\D'f\\'~o' 

A. 5. Interpolation Estimates. Here we record several interpolation estimate utilized in our proof. 
Since they have been proved in lemma A.6-A.8 and lemma 3.18 in [2], we will omit the proof now. 

Lemma A. 10. For s,q > and < r < s, we have the estimate 

(A.5.1) \\f\\HH^)^y\\^'-'^ikU^^^^^^^ 
whenever the right hand side is finite. 

Lemma A. 11. Let Vf be the Poisson integral of f , defined on S. Let A > 0, g, s e N, and r > 0. Then 
the following estimates hold. 

(1) Let 

(A.5.2) 6= 1-0= ^±A^ 

^ q+s+\ q+s+\ 

Then 

(A.5.3) l|V«P/f^o < (||/a/||'ho)'(||^)''+V||1o)' ' 

(2) Letr + s>l, 
(A.5.4) e= ' + '-\ 1-0= ^ + ^ + 1 



q + s + r + X q + s + r + \ 



Then 



(A.5.5) W^'VfWl^ < ||/a/||^o ) ( WD'^^'fWl. 

(3) Lets>l. Then 
(A.5.6) l|V«P/||i„o < WD'^fWl. 

Lemma A. 12. Let f be defined on S. Let A > 0. Then we have the following estimates. 
(1) Let q,s £ (0, oo) and 

(A.5.7) 0=^^ 1 



q + s + X q + s + X 

Then 



(A.5.8) IIDVIlLo < (|IWIIho)"(||^'+V||Lo' ' ' 
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(2) Letq,sGn,r>0,r + s>l, 

r + s-1 q + X + 1 



q + s + X q + s + A 



r+s-1 ^ ^ g+A+1 



q+s+r+X q+s+r+X 



(A.5.9) e X - " 

^ q+s+r+X q+s+r+X 

Then 

(A.5.10) WDOfWU <(^\\hff„o^'(^\\D<'+^f\\l.y ' 

Lemma A. 13. Let f be defined in Cl. Let X > 0, q,s G and r > 0. Then we have the following 

estim,ates. 

(1) Let 
(A.5.11) 

Then 

(A.5.12) IIDVIIho < (|Mho)'(||^'+V||^o)' 

(2) Letr + s>l 
(A.5.13) 

Then 

(A.5.14) IIDVIlioo < (l|/A/||^.)'(||i5«+V||?,.+.)' ' 

and 

(A.5.15) II^VIIloo(e) < ' 

A. 6. Continuity and Temporal Derivative. In the following, we give two important lemmas to 

connect L^H^ norm and L°^H^ norm. 

Lemma A.14. suppose that u e L2([o^ T]; if^i (f^)) and dfU e L2([0, T]; ff^^ (17)) for si > S2 > and 
s = {si+ S2)/2. Then u G C"([0, T]; H"{V,)) and satisfies the estimate 

(A.6.1) Ihlli^.^, < \\um\%. + llulli.^., + 

where the L^H'^ norm and L°°H'^ norm are evaluated in [0,r]. 

Proof. Considering the extension theorem in Sobolev space, we only need to prove this result in the i?" 
case. The periodic case can be derived in a similar fashion. Using Fourier transform, 

= 2/ {o''m,tMr\dtu{um< [ {o^''m,t)fd^+ [ {o'''\dtu{u)fd^ 

jRn jRn jRn 

= \\u{t)\\%s,+\\dtu{t)rHs. 

So integrate with respect to time on [0, t] 

+ \\dtu\\l 

□ 

The following lemma shows the estimate in another direction. 

Lemma A.15. For any u e Z/°°([0, T]; iJ'^) within [0,r], we must have u € L^([0, T]; ff'^) and satisfies 
the estimate 

(A.6.2) < TMl^H" 

Proof. The result is simply based on the definition of these two norms. □ 
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A. 7. Extension Theorem. The fohowing arc two extension theorems wliicli will be used to construct 
start point of iteration from initial data in proving wellposedness of Naiver-Stokes-transport system. 
Since it is identical as lemma A. 5 and A. 6 in [1], we omit the proof here. 

Lemma A. 16. Suppose that dlu{0) e H'^'^~'^^{fl) for j = 0,...,N, then there exists a extension u 
achieving the initial data, such that 

diu e L'^{[0, oo); H^^-^^+\fl)) n L°°{[0, oo); H^^-^^{n)) 

for j = 0, . . . ,N . Moreover, 

iV n n Af n 



(A.7.1) 



3=0 



^2^2JV-2j + l 



< 



j=0 



H2N-2j 



Lemma A. 17. Suppose that dlp{Q) G iJ^^ ^-^ ^(f}) for j = 0, . . . , N — 1, then there exists a extension 

p achieving the initial data, such that 

for j = 0,...,N— 1. Moreover, 

N-l .. o „ iV-l 



(A.7.2) 



j=0 



£,2JJ2N-23 



+ 



dip 



j=0 



U2N-2j-l 



Appendix B. Estimates for Fundamental Equations 

B.l. Transport Estimates. Let E be either infinite or periodic. Consider the equation 

. . ( dtV + u-Dr] = g in E x (0,T) 

\v{t = 0)=rjo 

We have the following estimate of the regularity solution to this equation, which is a particular case of a 
more general result proved in proposition 2.1 of [5]. Note that the result in [5] is stated for Y, = R^, but 
the same result holds in the periodic case as described in [9]. 

Lemma B.l. Let rj he a solution to equation (B.1.1). Then there exists a universal constant C > 
such that for any s > 3 

(B.1.2) <exp(^C^ ll^w(r)||^3/2dr^ ^ll^llna +^ IIsWIIh^ dr^ 

Proof. Use p = P2 = 2, N = 2 and S = s in proposition 2.1 of [9], along with the embedding H^/'^ M- 

Lemma B.2. Let r] be a solution to (B.1.1). Then there is a universal constant C > such that for 

any < s < 2 

(B.1.3) 

sup ||??(r)||jjs(s) < exp 



□ 



(B.2.1) 



Proof. The same as lemma A. 9 in [2]. 

B.2. Elliptic Estimates. We need two different type of elliptic estimates in our proof either in weak 
sense or strong sense. 

Lemma B.3. (weak solution) We call {u,p) G Hl{^l) x i?°(ri) a weak solution to the elliptic equation 

' -Au + Vp = (t)eH-^{n) in n 

V ■u = i} eH^{fl) in O 

u = on T, 

u — on Tib 

if'V-u = ip in n and {u,p) satisfy the weak formulation that for arbitrary f € -ffg (Q), it holds that 
(B.2.2) {m : D/)ffo + (p, V • /)ffo = (</., 

where (•, ■)ho represents the inner product in H^{n) and (•, denotes the dual pairing between Hq{Q) 
and H~^{fl). Then we have the estimate 

(B.2.3) \\u\\i. + wpwIo < ml-. + \mlo 
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Proof. An obvious modification for the proof of lemma 3.3 in [4] impHcs that for any tp G H'^{d), there 
exists awe H^^i}) such that V ■ v = ip in il, and satisfies the estimate \\v\\^i < \\4'\\'ho- Then we may 
change to the unknown w = u — v, which satisfies the system 



(B.2.4) 



—Aw + Vp = <p + Av in O 

V ■ w = in ft 

w = on T, 

w = on E(, 



Then we have the new weak formulation that for arbitrary / e Hq{CI), it is valid that 
(B.2.5) {Bw : ©/)ho + (p, V • /)ho = {<f>, - pv : ©/)ho 

We may further restrict that / satisfies V • / = as the pressureless weak formulation that 
(B.2.6) {Bw : ©/)ho = (</>, - {Bv : P/)ho 

Ricsz theorem implies that there exists a pressureless weak solution w G Hq (ft) satisfying V • w = and 

the estimate 

(B.2.7) \\w\\% < mi^, + Ml, < mi^, + M% 

Then it is well-known that pressure can be taken as the Lagrangian multiplier for the Navier-Stokes 
system. Similar to proposition 2.9 in [1], we have that there exists a. p G H^{Cl) such that the weak 
formulation (B.2.5) holds and it satisfies 

(B.2.8) ibii^o < Ml-, + \mlo + Ml, < \mi., + Mio 

Therefore, {u,p) satisfy the weak formulation (B.2.2) and our result easily follows. □ 
Lemma B.4. (strong solution) Suppose {u,p) solve the equation 

-Au + Vp = <j>G H'-'^ifl) in n 



(B.2.9) 



V ■u = ^ G H'-^n) in n 

(p7-Du)e3 = <^eiI'-3/2(S) on S 
u = on S(, 



in the strong sense. Then for r>2, 

(B.2.10) Ml. + Ml.., < \ml... + Ml.., + Mlr-s;. 

whenever the right hand side is finite. 

Proof. The same as lemma A. 15 in [1]. □ 
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